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PREFACE 


The  research  reported  herein  was  conducted  during  1981  and  1982  in 
the  Hydraulics  Laboratory  (HL)  of  the  US  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES).  The  investigations  were  funded  under  the  Environ¬ 
mental  and  Water  Quality  Operational  Studies  (EWQOS)  Program,  Task  IA.3 
"Improve  and  Verify  Physical  Hydrodynamic  Modeling  Techniques  for  Reser¬ 
voir,"  which  was  sponsored  by  the  Office,  Chief  of  Engineers  (OCE).  The 
OCE  Technical  Monitors  for  EWQOS  were  Dr.  John  Bushman  and  Messrs.  Earl 
Eiker  and  James  L.  Gottesman.  The  EWQOS  Program  was  assigned  to  the  WES 
under  the  direction  of  the  Environmental  Laboratory  (EL). 

Mr.  Henry  Simmons,  Chief  of  the  HL,  and  Mr.  John  L.  Grace,  Jr., 
Chief  of  the  Hydraulic  Structures  Division  (HS) ,  directed  the  research. 
Dr.  Dennis  R.  Smith,  former  Chief  of  the  Reservoir  Water  Quality 
Branch  (RWQB)  in  HS,  supervised  the  effort.  Messrs.  Mark  S.  Dortch 
and  Steven  C.  Wilhelms  conducted  the  studies  reported  herein  and  pre¬ 
pared  the  bulk  of  this  report.  Assisting  in  laboratory  work  were 
Messrs.  Calvin  Buie,  Alvin  Myers,  and  James  Daub.  Mr.  Jeffery  P. 
Holland,  Chief  of  the  RWQB,  reviewed  the  report  and  assisted  in  its 
final  preparation.  The  report  was  edited  by  Ms.  Jamie  W.  Leach  of 
the  WES  Publications  and  Graphic  Arts  Division. 

Dr.  Jerome  Mahloch  was  Program  Manager  of  the  EWQOS  during  these 
studies.  Dr.  John  Harrison  was  Chief  of  EL.  Director  of  WES  was 
COL  Allen  F.  Grum,  USA.  Technical  Director  was  Dr.  Robert  W.  Whalin. 

This  report  should  be  cited  as  follows: 

Dortch,  M.  S. ,  Wilhelms,  S.  C. ,  and  Holland,  J.  P.  1985. 

"Physical  Modeling  of  Reservoir  Hydrodynamics,"  Technical 

Report  E-85-14,  US  Army  Engineer  Waterways  Experiment 

Station,  Vicksburg,  Miss. 
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PHYSICAL  MODELING  OF  RESERVOIR  HYDRODYNAMICS 


PART  I:  INTRODUCTION 


Background 


1.  In  concern  for  environmental  quality,  water  resource  engineers 
and  planners  evaluate  water  quality  within  and  released  from  man-made 
lakes  or  reservoirs.  Typically  their  concerns  extend  over  the  distribu¬ 
tion  of  water  quality  characteristics  over  the  depth  of  the  reservoir 
and  the  quality  of  water  released  through  the  reservoir  flow  control 
system.  Some  of  the  quality  constituents  of  Interest  are  temperature, 
suspended  and  dissolved  solids,  dissolved  oxygen,  reduced  metals,  and 
other  water  quality  parameters.  Of  significant  importance  to  the  qual¬ 
ity  of  water  in  a  reservoir  is  the  transport  of  these  constituents  into, 
through,  and  out  of  the  reservoir. 

2.  A  variety  of  techniques  and  tools  to  evaluate  reservoir  water 
quality  are  available  to  the  water  resource  planner  or  engineer.  Many 
of  these  tools  Involve  transport  descriptions.  Physical  models  of 
varied  descriptions  can  be  used  to  evaluate  the  hydrodynamics  and  trans¬ 
port  phenomena  of  a  reservoir.  Physical  models  are  either  near-field 

or  far-field  models  and  are  used  to  evaluate  withdrawal  characteristics, 
inflow,  entrainment,  mixing,  dispersion,  dilution,  velocity  patterns, 
and  density  currents.  One-dimensional  numerical  water  quality  models 
are  used  extensively  to  evaluate  the  thermal  and  quality  regimes  of  a 
reservoir.  Two-dimensional  numerical  models  are  being  developed  for 
evaluating  hydrodynamics  and  thermal  and  quality  aspects  of  a  reservoir 
along  two  axes  (usually  the  vertical  and  longitudinal  axes)  simultane¬ 
ously.  Hybrid  modeling,  which  entails  a  combination  of  physical  and 
numerical  modeling,  can  be  used  to  provide  higher  quality  evaluations  of 
the  reservoir  environment. 

3.  Near-field  physical  models  reproduce  localized  hydrodynamics. 
Local  mixing,  entrainment,  and  withdrawal  flow  conditions  can  be  simu¬ 
lated  and  their  characteristics  can  be  defined.  For  hybrid  modeling. 
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mixing,  entrainment,  or  withdrawal  is  physically  simulated,  mathemati¬ 
cally  quantified,  and  incorporated  into  numerical  models.  The  near¬ 
field  physical  model  provides  a  site-specific  description  of  the  local¬ 
ized  hydrodynamics  that  result  in  a  more  accurate  description  of  boundary 
conditions  in  the  numerical  model. 

4.  Far-field  physical  models  simulate  the  overall  water  body  hy¬ 
drodynamics  extending  beyond  the  localized  flow  field  typically  investi¬ 
gated  in  a  near-field  model.  The  effects  of  complex  geomorphology  or 
operations  have  been  investigated  in  far-field  models.  Of  interest  in 
these  models  are  the  effects  of  local  mixing  upon  the  far-field  and  the 
velocity  patterns  and  movement  of  density  currents  through  the  reser¬ 
voir.  The  results  from  these  models  are  usually  incorporated  into  a 
numerical  model,  exemplifying  the  hybrid  modeling  approach. 

5.  Physical  and  numerical  models  typically  complement  each  other. 
The  information  gained  with  one  model  is  enhanced  by  the  knowledge  gained 
from  the  other.  However,  neither  numerical  nor  physical  models  can  be 
used  indiscriminately.  Numerical  models  are  no  more  accurate  than  the 
assumptions  associated  with  their  development  and  the  numerical  proce¬ 
dure  used  to  arrive  at  a  solution.  Physical  models  are  no  more  accurate 
than  their  scaling  criteria.  Some  physical  aspects  of  prototypes  are 
exaggerated  in  a  model;  for  instance,  bottom  roughness  may  be  dispropor¬ 
tionately  large  in  a  physical  model.  This  does  not  preclude  the  use  of 
physical  models.  Rather,  these  models  yield  excellent  results  if  geo¬ 
metric  and  hydrodynamic  similarity  is  preserved.  Additionally,  for  com¬ 
plex  boundary  or  operational  conditions,  physical  models  may  be  the  only 
practical  technique  to  describe  the  hydrodynamics. 

Purpose  and  Scope 

6.  The  objective  of  this  report  is  to  overview  the  applicability 
of  scaled  physical  models  for  the  solution  of  water  resources  problems. 
Particular  attention  will  be  given  to  the  evaluation/verification  of  the 
criteria  used  to  scale  various  reservoir  hydrodynamic  processes  in  phy¬ 
sical  models.  These  processes  include  selective  withdrawal,  density 
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currents,  hydraulic  jets,  and  mixing  processes.  Undistorted  and  dis¬ 
torted  physical  model  scaling  relationships  are  presented  with  emphasis 
on  criteria  for  model  distortion.  Example  applications  of  physical 
models  are  presented  and  distorted  model-prototype  comparisons  are  made. 
The  theoretical  basis  of  general  scaling  relationships  for  physical 
models,  although  well  established  and  documented  in  other  texts,  is  also 
presented  in  order  to  provide  the  reader  access  to  these  concepts  within 
a  single  document. 
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PART  II:  FROUDIAN  SCALING  FOR  PHYSICAL  MODELS 


7.  Dimensionless  groupings  that  characterize  stratified  flow  can 
be  deduced  from  the  Navier-Stokes  equations.  From  these  groupings, 
scaling  criteria  can  be  developed  for  use  in  the  application  of  physical 
modeling  to  water  resource  problems.  If  incompressibility  is  assumed, 
these  equations  may  be  expressed  in  vector  notation  as 


|^  +  q-Vq  =  g-  ^Vp  +  yV^q  (1) 

where 

q  =  ui  +  vj  +  wk,  a  three-dimensional  velocity  vector 
u,v,w  =  velocity  in  the  x- ,  y- ,  z-direction,  respectively 
t  =  time 

q  •  V  =  "dot"  product  of  the  velocity  vector  and  V 

Q  _  Q  _  0  _ 

V  =  gradient  operator, 

g  =  gravitational  acceleration  vector 
p,(j  =  fluid  density  and  viscosity,  respectively 

p  =  pressure  ?  ? 

2  3^  3^ 

V  =  Laplacian  operator,  — -  +  — -  +  —  2 

dx^  dy^  dx 

8.  Introducing  a  small  density  perturbation  p'  such  that 
p  =  p^  +  p'  (where  p^  is  the  density  of  a  homogeneous  fluid  and 
(p'/p^)  <<1)  will  result  in  approximating  incompressibility  and  will 
allow  the  effects  of  density  differences  to  appear  in  the  equation. 

Let  the  pressure  term  p  be  defined  .such  that 

P  =  +  P^g  •  r  +  p-  (2) 

where 

p  =  known  pressure  at  reference  location  c 
c 

r  =  vector  from  location  c  to  location  of  p 
p’  =  pressure  perturbation  from  homogeneous,  hydrostatic  fluid 


Substituting  for  p  ,  Equation  1  becomes 


at 


Vq  =:  g  - 


y  (p^ 


8  •  t  +  p  ’  ) 


P' 


V^q 


(3) 


The  coefficient  of  the  pressure  term  may  be  expanded  in  a  binomial 
series , 


Assuming  that  the  second-order  and  higher  order  terms  in  Equation  4  are 
approximately  zero  (which  is  reasonable  since  (p'/P^)  <<1)  and  since 

Vp^  =  0  (p^  by  definition  is  constant),  Vg  •  r  =  g  ,  and 

2 

(p'/P^)  Vp'  ~  0  (p'  and  p'  are  small),  the  pressure  term  (Equation  2) 
may  be  reduced  to 


8  +  ^  Vp'  -  ^  (5) 

^o  ^o 

Letting  p'  =  Ap  ,  some  small  incremental  density  difference,  and 
p/(p^  +  p')  =  V  ,  kinematic  viscosity,  and  substituting  Equation  5  into 
the  pressure  term  of  Equation  3,  the  e(|uation  of  fluid  motion  based  on 
the  Boussinesq  incompressibility  approximation  can  be  written  as 


3t 


q 


Vq  =  ^  8  -  ^  Vp'  +  vV^  q 


(6) 


9.  To  develop  similitude  criteria  for  modeling.  Equation  6  should 
be  nondimensionalized .  By  orienting  the  z-axis  in  the  vertical  direc¬ 
tion  (measured  positive  upward)  the  components  of  the  gravitational 
vector  g  become 


8  =  -8k 

where  k  is  the  unit  vector  in  the  z-direction.  The  z-component  of  Equa¬ 
tion  6  is 


i 


Equation  7  can  be  made  nondimensional  with  the  following  relationships: 


^  w  c 
w  =  u  .  £ 


X  =  L  .  y 


t  -S  U  -V 

L/U  ’  “  “  U  ’ 


Z  z  =  - 
L  ’  L 


where  the  " 
stants,  and 


denotes  a  dimensionless  variable,  p 


,  g  are  con- 


U  =  characteristic  velocity 
L  =  characteristic  length 

Substituting  the  above  expressions  into  Equation  7  and  simplifying  the 

2 

result  in  an  equation  that,  when  multiplied  by  L/U  ,  is  a  dimensionless 
mathematical  description  of  fluid  motion  in  the  z-direction: 


Inertial  Force 


Gravita-  Pres- 
tional  sure 
force  force 


Viscous 

Force 


A  similar  equation  may  be  written  for  the  x-  and  y-directions .  This 
equation  contains  two  dimensionless  parameters  that  are  familiar  to 
hydraul icians  and  hydrodynamicists :  the  Reynolds  number  (Re),  which  is 
the  inverse  of  the  coefficient  in  the  viscous  force  term;  and  the  densi 
metric  Froude  number  (Fr;  hereafter  referred  to  as  the  Froude  number), 
which  is  the  inverse  square  root  of  the  coefficient  in  the  gravitations 
force  term. 

10.  For  dynamic  similitude,  the  ratios  of  prototype  forces  to 
homologous  model  forces  must  be  equal.  That  is,  from  Equation  8,  the 
ratio  of  prototype  inertial  forces  to  model  inertial  forces  must  equal 


the  ratio  of  prototype-to-model  gravitational  forces,  which  must  equal 
the  ratio  of  prototype-to-model  viscous  forces,  which  must  equal  the 
ratio  of  prototype-to-model  pressure  forces.  Examination  of  Equation  8 
shows  that  the  equating  of  the  respective  prototype  and  model  Froude  and 
Reynolds  numbers  results  in  the  ratios  of  the  gravitational  and  viscous 
prototype-to-model  forces  being  equal  to  unity.  Since  the  equation  is 
dimensionless,  this  further  dictates  that  the  ratio  of  the  respective  in¬ 
ertial  forces  be  unity  as  well.  From  the  previous  definition  of  dynamic 
similitude,  the  ratio  of  prototype-to-model  pressure  forces  must  also  be 
unity.  Thus,  the  equating  of  respective  prototype  and  model  Froude  and 
Reynolds  numbers  ensures  the  existence  of  complete  dynamic  similitude. 

11.  It  is  usually  impossible  to  provide  complete  dynamic  simili¬ 
tude  since  scale  reductions  affect  some  dimensionless  ratios.  However, 
one  or  more  of  the  forces  may  be  absent  or  negligible  for  a  given  flow 
regime.  In  fluid  systems  with  a  free  surface,  such  as  open  hydraulic 
channels,  and  water  bodies,  such  as  reservoirs,  inertial,  pressure,  and 
gravitational  forces  often  dominate.  Thus,  the  Froude  number  should  be 
preserved  when  scaling  from  prototype  to  model. 

12.  However,  although  gravitational  forces  are  of  major  impor¬ 
tance  when  modeling  open  channel  hydraulics,  viscous  and  surface  tension 
forces  are  not  always  negligible.  For  reservoirs  and  other  large  water 
bodies,  viscous  forces  may  be  important  in  some  flow  regimes.  Since 
water  is  used  in  most  models  of  hydraulic  systems,  it  is  not  possible  to 
simultaneously  preserve  Froude  and  Reynolds  number  equalities.  Experi¬ 
ence  has  shown  that  if  the  fundamental  character  of  flow  is  preserved  in 
both  prototype  and  model  (i.e.,  both  systems  are  turbulent),  then  simil¬ 
itude  relationships,  valid  for  transfer  of  quantities  from  model  to  pro¬ 
totype,  can  be  developed  based  on  the  Froude  number  only  for  open  channel 
hydraulics.  Limitations  and  requirements  for  Freudian  modeling  when  vis¬ 
cous  forces  are  important  are  discussed  in  Part  III. 

Scaling  Relationships 

13.  The  Freudian  scaling  criterion  for  density  stratified  flow 


.y 

*  H  » 


-V  . 


dictates  that  the  Fronde  number  be  equal  in  model  and  prototype.  That  is 


or 


Fr  =  Fr 
m  p 


(9) 


where 


m,p 

V  ,V 
m  p 

Ap  ,  Ap 
^m  p 

P  .  P 

m  '^p 


L  ,  L 
m  p 


=  model  and  prototype,  respectively 
=  velocity  in  model  and  prototype,  respectively 
=  density  difference,  model  and  prototype,  respectively 
=  density  in  model  and  prototype,  respectively 
=  gravitational  acceleration 

=  characteristic  length  in  model  and  prototype, 
respectively 


By  providing  identical  density  differences  in  model  and  prototype  such 
that 


_ m  _  p 

P  ~  P 

m  p 

the  Freudian  scaling  relations  for  velocity,  time,  and  discharge  can 
be  simplified.  This  relationship  is  assumed  in  subsequent  modeling 
criteria . 

14.  Undistorted  near-field  reservoir  models  have  equal  horizontal 
and  vertical  scales.  These  models  are  geometrically  and  hydrodynami- 
cally  similar  to  the  prototype.  Distorted  scale  reservoir  models  typi¬ 
cally  have  a  vertical  scale  that  is  exaggerated  with  respect  to  the  hor¬ 
izontal  scale  to  allow  sufficient  depth  to  produce  turbulent  flow  and 
permit  stratified  flow  visualization.  This  results  in  different  scaling 
relationships  for  undistorted  and  distorted  models. 


Undistorted  modelin; 


15.  Several  scaling  relationships  are  required  for  the  design  of 
an  undistorted  model.  These  include  length,  velocity,  time,  and  dis- 
char^.  The  velocity  scale  relation  is  determined  by  rearranging  Equa¬ 
tion  9  into  a  velocity  ratio  (denoted  by  the  upper  case  R  subscript) 
for  model  and  prototype. 


V 


R 


(10a) 


Cancelling  terms  results  in 


(10b) 


The  velocity  scaling  ratio  is  the  square  root  of  the  length  ratio.  The 
time  scale  relationship  is  determined  by  substituting  velocity  L/t  for 
V  in  Equation  9.  Cancelling  terms  again  results  in  a  time  scale  ratio 


t 


R 


(11) 


where  t^,  t^  =  time  in  the  model  and  prototype,  respectively.  To  deter¬ 
mine  the  discharge  scaling  function,  the  continuity  relationship  is  used: 

Q  =  AV 

where 

Q  =  volumetric  flow  rate 
A  =  cross-sectional  area  of  flow 
V  =  velocity  of  flow 

Solving  for  velocity  from  the  continuity  equation  and  substituting  into 
Equation  9  results  in  the  discharge  scaling  relationship: 


Distorted  modeling 

16.  Scale  relations  for  distorted  scale  physical  models  can  be  de 
veloped  in  an  analogous  manner  with  slight  differences.  The  horizontal 
and  vertical  length  scale  ratios  for  a  distorted  model  are  defined  as 


(13a) 


Y„  =  Y  /Y 
R  pm 


(13b) 


where 


Xjj,  Yj^  =  length  scale  ratios  in  horizontal  and  vertical  directions 
respectively 

X  ,  Y  =  prototype  length  in  horizontal  and  vertical  directions, 


P  P 


respectively 


X  ,  Y  =  length  of  model  in  horizontal  and  vertical  directions, 
m  m  •  1 

respectively 

In  a  distorted  model,  the  scaling  ratio  for  area  in  a  horizontal  plane 
A  Is 


A  =  x; 

Xr  r 


(14a) 


while  in  a  vertical  plane  the  area  scaling  ratio  is 


(14b) 


The  characteristic  length  in  the  gravitational  term  (denominator)  of 
Equation  9  is  in  the  y-direction,  the  direction  of  gravity.  Therefore, 


Dimension 


Ratio 


Undistorted  Model 

Length 

=  L 

P 

m 

Velocity 

=  ^R 

/2 

Time 

'^R 

=  "r 

/2 

Volume  flow  rate 

Qr 

/2 

Density  difference 

"^R 

=  1 

Distorted  Model 

Length  in  the  horizontal 

X 

m 

Length  in  the  vertical 

Yr  = 

’'p' 

Y 

m 

Area  in  the  horizontal  plane 

A 

^R 

II 

Area  in  the  vertical  plane 

s 

=  ^R 

"r 

Velocity 

II 

> 

yl/ 

R 

2 

Time 

"r  = 

Volume  flow  rate 

II 

06 

O' 

V 

1 

Density  difference 

=  1 

18.  Model-to-prototype  length  scale  ratios  for  an  undistorted 
model  may  range  from  1:5  to  1:200,  depending  on  the  area  and  flow  regime 
of  interest.  Distorted  scales  may  range  from  1:800  to  1:4,000  in  the 
horizontal  (X^^)  and  1:20  to  1:200  (Y^^)  in  the  vertical.  Typically, 
distortion  factors  i^^nge  from  4  to  40. 


15 


PART  III:  DISCUSSION  OF  THE  ROLE  OF  DISTORTION  IN 
RESERVOIR  MODELING 


Criteria  for  Distortion 


19.  When  modeling  a  large  open  water  body  such  as  a  surface  res¬ 
ervoir,  there  may  be  several  flow  regimes  of  interest  such  as  those 
dominated  by  inertia,  buoyancy,  and  turbulent  dispersion.  One  may  be 
interested  in  not  only  what  occurs  near  a  perturbation,  but  also  what 
occurs  at  some  distance  (far  field).  Modeling  for  several  interests 
simultaneously  may  require  conflicting  modeling  criteria.  Freudian 
laws  are  still  required  to  give  the  proper  balance  between  inertial  and 
gravitational  forces,  but  vertical  length  scale  exaggeration  or  distor¬ 
tion  may  be  necessary  to  provide  adequate  turbulence  and  proper  viscous 
effects . 

20.  Modeling  such  large  water  bodies  often  necessitates  for  econ 

omy  of  size  a  large  horizontal  length  scale  ratio  X_  where 

K 


m 

and  Xp  and  X^  are  the  horizontal  dimensions  of  the  prototype  and 
model,  respectively.  With  a  large  horizontal  scale  ratio,  a  smaller 
vertical  scale  ratio  Y  must  usually  be  selected  to  provide  adequate 
model  depth  for  data  collection  and  turbulent  flow.  Thus,  the  vertica 
scale  is  exaggerated  or  distorted  and  the  distortion  factor  DF  is 


1 


The  concerns  related  to  this  modeling  are  how  much  distortion  is  needed 
and  what  is  too  much.  The  answers  to  these  questions  depend  on  the 
physical  processes  under  study  and  an  understanding  of  the  modeling  fea 
tures  that  may  be  lost  or  compromised  for  a  given  factor  of  distortion. 

21.  Physical  processes  that  may  be  of  interest  in  a  reservoir 
model  include  inflow,  outflow,  density  currents,  turbulent  mixing  and 


dispersion,  buoyant  jets,  and  heat  exchange.  These  processes  are  af¬ 
fected  by  the  geometry,  densimetric  Froude  number,  Reynolds  number, 
bottom  roughness,  and  coefficient  of  surface  heat  exchange.  A  discus¬ 
sion  of  the  modeling  of  these  processes  is  presented  in  the  context  of 
distorted  reservoir  model  applicability  based  on  equality  of  the  densi¬ 
metric  Froude  number  and  density  differences  and  approximately  equal 
density  and  viscosity  (same  fluid,  water)  in  the  prototype  and  model 

(p  =  Ap  =  p  =  1.0  for  the  prototype/model  ratios  of  density,  density 
K  K  K 

difference,  and  viscosity,  respectively). 

Inflow  and  outflow 

22.  Inflow  and  outflow  are  boundary  conditions  for  the  resulting 
internal  currents  in  density  stratified  reservoirs.  The  volumetric  flow 
rates  for  inflow  and  outflow  are  scaled  according  to  Froudian  scaling 

as  presented  in  paragraph  16.  This  scaling  is  appropriate  for  inflow 
since  inertia  and  buoyancy  will  be  properly  accounted  for  in  the  model. 

23.  In  some  instances,  simply  using  distorted  Froudian  scaling 
criteria  may  not  result  in  simulating  the  vertical  extent  of  the  zone  of 
withdrawal  during  outflow.  Many  reservoir  outlets  can  be  well  described 
as  point  sinks  (Bohan  and  Grace  1973).  The  vertical  distance  from  the 
elevation  of  the  sink  to  the  upper  or  lower  withdrawal  limit  Z  is 
dependent  upon  the  discharge  Q  and  the  density  difference  between  the 
elevations  of  the  sink  and  the  withdrawal  limit  Ap  This  distance 
can  be  approximated  by 


Z  oc 


(20) 


The  limits  of  withdrawal  should  scale  according  to  the  vertical  length 

scale  such  that  Z  =  Y  where  Z  equals  the  ratio  of  prototype  limit 
K  K  *^3/2 

to  model  limit.  Substituting  Q_  =  X„Y  (Froudian  discharge  scaling 

K  K  A 

from  Equation  19)  and  p„  =  Ap_  =1.0  into  the  above  equation,  it  is 

^  ^  2/5 

obvious  that  will  not  equal  Y  but  (OF)  Y„  .  Therefore,  for 

K  K  n 

the  general  case  of  a  long  reservoir  (DF  >  1),  the  model  withdrawal  zone 
would  not  extend  as  far  vertically  as  the  prototype  withdrawal  zone.  If 


the  outlet  is  a  rectangular  orifice,  then  direct  horizontal  and  vertical 


scaling  results  in  a  vertical  model  dimension  stretched  with  respect  to 
the  horizontal  dimension.  This  allows  some  remedy  to  the  dilemma  since 
such  stretched  outlets  may  cause  a  deepening  of  the  withdrawal  zone  be¬ 
yond  that  described  for  a  point  sink.  To  simulate  point  sink  with¬ 
drawal,  it  may  be  necessary  to  adjust  the  vertical  extent  of  the  outlet 
opening  to  simulate  the  withdrawal  limits  for  the  specified  scaled 
discharge . 

24.  When  a  prototype  release  structure  can  be  catagorized  as  dis¬ 
charge  through  a  line  sink,  adjustments  to  the  model  outlet  should  not 
be  necessary.  This  is  because 


where  q  is  the  unit  discharge.  Using  the  discharge  scaling  stated  in 
paragraph  23,  it  is  apparent  that  the  above  equation  results  in 
=  Yp  .  It  should  be  emphasized  that  distorted  reservoir  models  are  not 
intended  for  comprehensive  analysis  of  the  near-field  withdrawal  distri¬ 
bution.  This  typ  ;  of  analysis  is  usually  conducted  with  an  undistorted 


model  with  a  length  scale  ratio  L  equal  to  both  X_  and  Y  . 

K  K  K 


Density  currents 


25.  Even  if  the  inflow-outflow  boundary  conditions  have  been  pro¬ 
perly  established,  the  resulting  model  density  currents  (convective 
spreatling)  must  be  evaluated  for  similitude.  The  interfacial  shear 
stress  developed  by  a  layer  of  fluid  flowing  above  or  beneath  another 
layer  of  fluid  may  be  disproportionately  large  in  an  undistorted  model 


compared  with  inertial  and  gravitational  forces.  The  bottom  shear  may 


also  be  disproportionately  large  in  an  undistorted  model  if  the  model 
roughness  is  too  great.  Vertical  scale  exaggeration  in  a  distorted  model 
is  often  used  to  compensate  for  excess  shear  and  bottom  roughness. 

26.  Proper  adjustment  for  bottom  roughness  is  required  when  at¬ 
tempting  to  model  density  currents  flowing  along  the  bottom  (underflows) 
and  to  provide  similar  damping  of  internal  waves.  In  an  undistorted 
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model,  the  scaling  ratio  for  the  Manning's  n  roughness  factor  is 


"R  =  ’-r” 

If  the  prototype  bottom  is  fairly  smooth,  as  with  a  sandy  bottom,  it  may 
be  impossible  to  decrease  the  model  n  sufficiently  since  the  model  boun¬ 
dary  may  already  be  hydrodynamically  smooth.  However,  this  may  be  over¬ 
come  through  vertical  scale  exaggeration  where  the  scale  relation  for 
n  is 


yl/6 

R 

Vdf 


(23) 


Thus,  it  may  be  possible  to  provide  the  proper  roughness  effect  by  using 

a  smooth  model  boundary  and  adjusting  the  distortion  factor.  Usually  it 

1/3 

is  desirable  to  have  DF  >  Y  ;  therefore,  n  must  be  greater  than 

K  m 

n^  .  This,  in  turn,  requires  exaggerated  roughness  elements  and  cali¬ 
bration  procedures  for  proper  physical  model  adjustment. 

27.  When  the  bottom  is  relatively  flat,  bottom  roughness  adjust¬ 
ment  alone  may  be  sufficient  for  modeling  bottom  density  currents.  When 
modeling  river  or  channel  flow,  it  may  be  necessary  to  adjust  the  bed 
slope  and  roughness  in  a  distorted- length  scale  model  to  give  the  proper 
stage-discharge  relationship.  With  a  reservoir  model,  the  water  surface 
is  flat  and  the  depth  may  vary  considerably  from  the  headwater  to  the 
dam.  If  the  model  is  to  reproduce  the  capacity-elevation  characteristics 
of  the  prototype,  then  it  would  not  be  possible  to  reduce  the  bed  slope 
to  account  for  the  excessive  steepness  caused  by  distortion.  In  some 
instances  this  may  necessitate  an  adjustment  to  the  bed  roughness  to 
compensate  for  the  excessive  slope. 

28.  Interflow  and  overflow  (surface  flow)  density  currents  are 
affected  by  interfacial  shear.  A  discussion  of  interflow  and  overflow 
phenomena  can  be  found  in  Appendix  A.  Because  nonviscous  intrusion 
(interflow  and  overflow)  characteristics  are  dependent  on  I roude  number, 
Freudian  scaling  (distorted  or  undistorted)  should  properly  simulate  a 


nonviscous  prototype  intrusion  if  the  intrusion  in  the  model  is  also  in 
the  nonviscous  range.  However,  length  scale  distortion  may  be  required 
to  ensure  that  a  nonviscous  intrusion  will  exist  in  the  model.  To  il¬ 
lustrate  this,  it  is  first  assumed  that  the  density,  density  difference, 
and  viscosity  of  the  model  equal  that  of  the  prototype.  Then,  if  the 
criteria  from  Appendix  A  for  nonviscous  spreading  are  maintained  in  the 
model  when  there  is  nonviscous  spreading  in  the  prototype,  from 
Equation  A9 


1.0 


(24) 


where 

=  the  ratio  of  prototype  inflow  intrusion  length  to  analogous 
model  length 

=  the  ratio  of  total  prototype  inflow  thickness  to  analogous 
R  model  thickness 


which  is  the  required  relationship  to  maintain  properly  scaled  viscous 
forces.  An  undistorted  model  would  require  L  <  1.0  to  satisfy  the 
above  relation,  which  means  the  model  would  have  to  be  at  least  as  large 
as  the  prototype.  A  distorted  model  would  require 

Yr  1  (DF)^/^ 


9 

r 


The  distorted-scale  requirements  are  not  as  restrictive  as  the  undis¬ 
torted,  although  Y  may  still  be  prohibitive.  Rearranging  the  above 
relationship,  the  minimum  distortion  required  to  provide  complete  simil¬ 
itude  for  convective  spreading  is  dependent  on  the  vertical  or  horizon¬ 
tal  scale  selected,  such  that 

DF  >  (25a) 

K 


o  r 


DF  >  X 


O.b 

R 


(25b) 


29.  Providing  nonviscous  spreading  in  the  model  for  the  complete 

range  of  nonviscous  spreading  in  the  prototype  may  require  fairly  large 

distortion  factors.  For  example,  if  a  horizontal  length  scale  ratio  of 

1,000  is  selected.  Equation  25b  calls  for  DF  >  63  which  results  in  Y 

~  K 

=  15.85  .  This  amount  of  distortion  and  the  resulting  Y_  may  be  unde- 

K 

sirable  and  unnecessary.  A  larger  vertical  length  scale  ratio  Y 

K 

may  be  more  practical  from  a  model  construction  and  operation  perspec¬ 
tive.  Additionally,  if  too  much  distortion  is  used,  wall  and  end  ef¬ 
fects  may  be  introduced.  Further,  excessive  distortion  can  result  in 
excessive  horizontal  (longitudinal  and  lateral)  dispersion  (which  will 
be  discussed  later).  To  avoid  such  problems,  a  distortion  factor  less 
than  that  specified  by  Equation  25b  is  accepted,  often  without  compro¬ 
mising  similitude  for  the  nonviscous  density  currents  of  interest.  For 
example,  the  nonviscous  spreading  zone  may  extend  beyond  the  area  of 
interest.  The  distortion  should  be  sufficiently  large  to  ensure  that 
the  prototype  nonviscous  spreading  distance  of  interest  remains 

nonviscous  for  the  scaled  distance  S,  in  the  model.  Existence  of  the 

m 

nonviscous  condition  can  be  estimated  for  the  prototype  and  model  for 
interflow  and  overflow  conditions.  However,  a  viscous  intrusion  slows 
down  and  becomes  thicker  with  time  as  buoyant  and  viscous  forces  domi¬ 
nate.  Froudian  modeling  criteria  would  therefore  no  longer  be  applica¬ 
ble  for  similitude.  Thus,  similitude  for  a  Froudian  model  may  lapse 
when  viscous  convective  spreading  occurs.  These  points  are  discussed  in 
greater  detail  in  Appendix  A. 

Far-field  dispersion 

30.  As  indicated  earlier,  Froudian  scaling  will  not  preserve  the 
prototype-model  Reynolds  number.  The  undistorted  model  Reynolds  number 
will  be  less  than  the  prototype.  However,  for  many  near-field  flow  prob 
lems,  the  flow  is  very  turbulent  and  the  area  of  interest  is  of  short 
length  (such  as  a  hydraulic  jump).  For  these  cases,  the  dissipation 

of  energy  is  mostly  due  to  turbulent  fluctuations  rather  than  laminar 
viscous  effects.  Turbulent  fluctuations  are  proportional  to  the  square 
of  the  average  velocity,  as  are  the  inertial  forces;  thus,  dissipative 
forces  associated  with  short,  very  turbulent  flows  are  similar  in 


undistorted  Froiidian  models.  Therefore,  near-field  turbulent  diffusion, 
which  is  also  related  to  turbulent  fluctuations,  can  be  simulated  with 
undistorted  Freudian  models.  For  far-field  dispersion  problems,  the 
flow  patterns  most  likely  will  not  be  short  and  not  as  turbulent  as  near¬ 
field  problems.  This  makes  it  more  difficult  to  account  for  dissipative 
forces  and  tu'bulent  fluxes.  Viscous  forces  and  laminar  diffusion  in 
the  model  may  dominate  the  turbulent  fluctuations.  With  an  uridistorted 
far-field  model,  a  reasonable  length  scale  would  most  likely  result  in 
such  a  shallow  depth  that  turbulent  flow  may  not  occur  in  the  model. 

Thus,  the  vertical  scale  is  exaggerated  and  a  distorted  model  is  used  to 
attain  turbulent  flow.  From  Froudian  law  with  length  scale  distortion, 
the  required  scale  ratios  for  proper  horizontal  (transverse)  and  verti¬ 
cal  dispersion  are 
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(26a) 


(26b) 


31.  The  actual  transverse  dispersion  coefficient  Dx  for  an  open 
channel  is  proportional  to  the  product  of  the  shear  velocity  u*  and 
the  depth  d  .  For  a  wide  channel 


oc  diu  =  d  -\/gds^ 


where  s  is  the  energy  slope.  Therefore,  the  prototype-to-model  ratio 
for  transverse  dispersion  using  Kquation  27  is 


‘r  Vdf 


Thus  it  is  apparent  that  the  transverse  model  dispersion  may  not  be 
scaled  appropriately  by  Equation  28.  Consider  that  density  stratifica¬ 
tion  and  turbulence  affect  the  vertical  dispersion.  These  effects  have 
been  represented  as 

D  =  D  f(R.)  (29) 

y  y^ 


where 


D  =  diffusion  coefficient  related  to  velocity  and  depth  in  a 
^o  nonstratif ied  reservoir 


f(R.)  =  a  function  of  Richardson  number  R. 
1  1 


Freudian  scaling  criteria  preserve  the  prototype/model  (R^  is  ap¬ 


proximately  equal  to  the  inverse  of  the  Froude  number  squared).  In  gen¬ 


eral,  D  a  ud,  where  u  is  horizontal  velocity.  Thus  the  vertical 


scale  ratio  for  dispersion  takes  the  form 

.1.5 
'R 


D  ay; 
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(30) 


Again  model  dispersion  scales  differently  from  that  desired  (Equa¬ 
tion  26b) . 

32.  In  the  distorted  model,  dispersion  can  be  physically  related 
to  the  scale  ratios  by  Equations  28  and  30.  The  desired  scale  relations 
for  dispersion  in  the  distorted  model  are  defined  by  Equations  26a  and 
26b.  The  desired  and  actual  distorted  scale  relationships  for  disper¬ 
sion  are  summarized  below: 


Desired  Distorted  Scale 
Direction  Ratio  for  Dispersion 


Horizontal 

(Transverse) 


Vertical 


DF 

K 


DF 


Actual  Distorted  Scale 
Ratio  for  Dispersion 


D 

^R 
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From  examination  of  the  above  relationships,  it  is  evident  that  turbu¬ 
lent  dispersion  in  the  distorted  model  is  too  fast  (overpredicted)  in 
the  horizontal  (transverse)  directions  and  too  slow  (underpredicted)  in 
the  vertical  direction.  Evidence  of  this  result  was  reported  by  Crick- 
more  (1972)  (as  referenced  on  p  304  in  Mixing  in  Inland  and  Coastal 
Waters ;  Fischer  et  al.  1979)  with  a  comparison  of  model  and  prototype 
observations  of  dispersion  of  a  continuous  source  tracer.  The  model 
plume  spread  faster  horizontally  and  slower  vertically  than  the  proto¬ 
type  plume.  Similar  results  for  lateral  dispersion  were  found  by  Moretti 
and  McLaughlin  (1974)  during  laboratory  reservoir  inflow  experiments  de¬ 
signed  to  determine  the  effects  of  distortion.  With  careful  adjustment 
of  model  rougfmess  (flow  inhibitors),  it  is  possible  to  alter  transverse 
and  vertical  mixing  to  improve  far-field  dispersion  modeling.  Despite 
possible  inaccuracies  in  dispersion  modeling,  distortion  is  required 
when  modeling  large  water  bodies  in  order  to  preserve  turbulent  flow  and 
to  reduce  shear  stresses  for  nonviscous  spreading  of  density  currents. 
Buoyant  jets/plumes 

33.  Source  momentum,  buoyancy,  and  geometry  are  all  important 
when  modeling  the  characteristics  of  a  free  turbulent  shear  flow  jet/ 
plume  (here  simply  jet  unless  otherwise  noted).  The  turbulent  shear 
flow  causes  the  entrainment  of  ambient  fluid,  thereby  increasing  the 
volume  flux  of  the  jet  along  its  trajectory.  If  an  undistorted  model 
jet  is  sufficiently  turbulent,  similitude  can  be  preserved  through 
Froudian  scale  relationships  (see  paragraph  30).  However,  a  distorted 
model  alters  the  model  source  geometry  from  that  of  the  prototype.  This 
alteration  may  change  other  jet  characteristics  such  as  entrainment. 

Thus,  jet  modeling  should  be  accomplished  without  distortion  in  a  near¬ 
field  model.  How  can  the  proper  effects  of  an  undistorted  jet  model  be 
created  in  a  distorted  scale  model? 

34.  If  a  line  diffuser  (line  source)  is  being  modeled,  then 
methods  outlined  in  Chapter  10  of  Mixing  in  Inland  and  Coastal  Waters 
(Fischer  et  al.  1979)  can  be  used  to  size  the  diffuser  in  the  distorted 
moilel.  The  diffuser  port  sizes  and  spacing  are  kept  undistorted  by  scal¬ 
ing  according  to  the  vertical  scale,  and  the  number  of  model  ports  is 
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reduced  by  the  distortion  factor.  Scaling  procedures  are  not  as  simple 
when  a  point  or  single  source  diffuser  and  jet  is  being  modeled.  The 
jet  entrainment,  dilution,  and  volume  flux  characteristics  must  be  known 
(possibly  from  results  of  a  near-field  undistorted  model),  and  the  dis¬ 
torted  model  source  geometry  must  be  adjusted  to  yield  the  same  charac¬ 
teristics  in  the  distorted  model.  If  the  jet  behavior  can  be  estimated 
from  relationships  developed  from  previous  analytical  and  experimental 
general  investigations,  then  it  may  be  possible  to  design  a  distorted 
model  jet  that  will  give  similar  behavior.  An  example  of  such  a  proce¬ 
dure  is  presented  in  Appendix  B. 

Heat  exchange 

35.  Reservoir  and  lake  water  temperatures  are  primarily  driven 

by  hydrological  and  meteorological  conditions  and  display  seasonal  vari¬ 
ations  (unless  the  lake  is  used  for  cooling  thermal  discharges).  It 
would  be  impractical,  if  not  impossible,  to  model  heat  exchange  char¬ 
acteristics  over  a  season  in  a  physical  model.  If  the  model  simulation 
time  is  not  too  long  (several  weeks),  meteorological  effects  (such  as 
solar  heating)  can  be  neglected,  and  the  model  can  be  used  to  analyze 
the  reservoir  hydrodynamics  for  the  event  of  interest. 

36.  Modeling  receiving  waters  for  thermal  discharges  (as  from  a 
power  plant)  may  require  modeling  surface  heat  exchange  if  the  receiving 
water  of  interest  covers  a  very  large  area.  Modeling  surface  heat  ex¬ 
change  requires 


=  1.0 


where 

k_  =  surface  heat  exchange  coefficient  ratio 
K  2 

K  =  surface  heat  exchange  coefficient,  Btu/ft  /“F/sec 

3 

y  =  specific  weight  of  water,  Ib/ft 
Cp  =  specific  heat  of  water,  Btu/lb/°F 
u  =  characteristic  velocity,  ft/sec 
With  =  1.0  ,  similitude  requires  Kj^  .  Thus, 

be  less  than  K  by  the  ratio  -v/L  .  This  can  be  achieved  by 
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controlling  the  atmosphere  in  the  laboratory.  If  the  laboratory  does 
not  have  atmospheric  control,  it  may  be  necessary  to  use  distortion  to 
compensate  for  the  excessive  heat  loss  rate  in  the  model.  Methods  for 
determining  the  amount  of  distortion  needed  are  discussed  in  Chapter  10 
of  Mixing  in  Inland  and  Coastal  Waters  (Fischer  et  al.  1979). 

Model-Prototype  Comparison 


37.  Tests  were  conducted  at  WES  to  compare  prototype  and  distorted 
model  behavior.  Two  lake  hydrodynamic  processes  were  studied.  First, 
interflow  and  overflow  (surface  flow)  density  currents  were  examined  with 
"prototypes"  and  models.  Then,  hydraulic  (water  jet)  destratification 
was  tested  for  "prototype"  and  model  conditions. 

Density  currents 

38.  Two  studies  were  conducted  to  compare  distorted  model  versus 
"prototype"  density  current  observations.  One  model  comparison  was  for 
DeGray  Lake,  which  is  a  Corps  lake  in  Arkansas.  The  other  comparison  was 
for  a  distorted  model  of  the  80.0-ft-long*  general  reservoir  flume  (re¬ 
ferred  to  as  the  General  Reservoir  Hydrodynamic  (GRH)  Facility)  located 
in  the  Hydraulics  Laboratory  of  WES.  In  this  case,  the  GRH  was  the 
"prototype . " 

39.  DeGray  Lake  model  study.  Flood  events  on  the  Caddo  River 
were  simulated  in  a  distorted  scale  physical  model  of  DeGray  Lake.  The 
flood  events  had  been  observed  (Ford,  Johnson,  and  Monusmith  1980)  and 
tracked  in  the  prototype.  Model  and  prototype  observations  were  compared 
to  identify  similarities  and  disparities  and  areas  of  modeling  theory, 
techniques,  and  procedures  that  should  be  improved. 

40.  The  model  had  a  horizontal  scale  ratio  of  2,400  and  a  verti¬ 
cal  scale  ratio  of  120,  with  a  distortion  factor  o-'  20.  Based  on  Freud¬ 
ian  scaling  criteria  (distorted  model  scaling  relationships  in  para¬ 
graph  17)  the  following  scale  ratios  may  be  used  to  quantitatively  trans¬ 
fer  model  (lata  to  prototype  equivalents: 

A  table  of  factors  for  converting  non-SI  units  of  measurment  to  SI 
(metric)  is  presented  on  page  3. 


Length  in  the  horizontal 

2,400 

Length  in  the  vertical 

120 

Area  in  the  horizontal 

5,760,000 

Area  in  the  vertical 

288,000 

Velocity 

10.9545 

Time 

219.089 

Volume  flow  rate 

3,154,881.934 

Density 

1 

The  model  was  constructed  of  clear  plastic  to  permit  flow  visualization. 
The  elevation-storage  relationship  and  general  lake  geometry  were  repro¬ 
duced.  No  attempt  was  made  to  simulate  roughness  of  the  lake  bottom  or 
sides . 

41.  Three  flood  events  were  studied  in  the  model:  29  March  1979, 
14  November  1978,  and  24  October  1976.  Plate  1  shows  the  density  struc¬ 
ture  of  the  lake  for  each  of  the  three  events.  The  density  profiles 
were  based  on  the  thermal  structure  of  the  lake  at  the  time  of  each  flood 
event.  These  observed  density  differences  between  the  lake  epilimnion 
and  hypolimnion  and  epilimnion  and  inflow  were  simulated  with  fresh  and 
saline  waters.  Plate  1  also  shows  the  stratification  patterns  used  in 
the  physical  simulations.  As  indicated  in  Plate  1,  the  model  stratifi¬ 
cation,  although  not  exact  replicates,  did  approximate  prototype  strati¬ 
fication.  Recreating  density  differences  is  relatively  easy,  but  repro¬ 
ducing  the  exact  prototype  density  distribution  in  the  model  is  much 
more  difficult.  Plates  2-4  show  the  flood  hydrographs.  The  hydrographs 
were  reproduced  by  controlling  inflow  with  a  manually  operated  valve. 

The  movement  of  the  inflow  was  tracked  in  the  model  by  dying  the  inflow 
and  marking  the  leading  edge  at  specific  times  during  each  test.  The 
inflow  movement  was  tracked  in  the  prototype  (Ford,  Johnson,  and  Monu- 
smith  1980)  by  mapping  the  movement  of  inflow  turbidity,  temperature, 
and  other  naturally  occurring  tracers  of  fluorescent  dye  that  was 


injected  into  the  river  upstream  of  the  lake. 

42.  Model-to-prototype  comparisons  were  made  by  comparing  travel 
time  and  distance  measured  along  the  mainstem  of  the  Caddo  River.  To 
ensure  a  realistic  comparison  of  model  and  prototype  travel  times,  time 
zero  was  selected  to  be  the  earliest  prototype  data  point  location. 
Subsequent  prototype  observations  were  referenced  to  this  time  zero. 
Similarly,  model  data  (after  conversion  to  prototype  equivalents)  were 
also  referenced  to  this  time-zero  location.  This  negated  any  model- 
prototype  differences  in  inflow  timing. 

43.  Figure  1  shows  results  of  tests  of  the  29  March  1979  flood 
event.  The  model  reproduced  the  movement  of  the  density  current  for 
this  event  very  well.  Equation  A6  was  evaluated  with  the  following  pro¬ 
totype  data:  h  ^  33  ft  ,  Ap  ~  0.00105  g/cc  ,  p  =  1.0000  g/cc  ,  v 
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=  10  ft  /sec  ,  and  £  2;  24,300  ft  .  This  yielded  a  value  of  473  which 
indicated  that  the  prototype  is  in  the  nonviscous  flow  regime  (see 
Appendix  A).  Model  data  of  h.^.  ~  0.3  ft  ,  Ap.j,  =  0.00085  g/cc  ,  p 
=  1.0000  g/cc  ,  V  =  10  ^  ft^/sec  ,  and  £  2  8.3  ft  yielded  a  value  of 
about  98,  indicating  that  model  flow  is  also  in  the  nonviscous  flow  re¬ 
gime.  Since  the  model  and  prototype  were  in  the  nonviscous  flow  regime, 
the  model  reproduced  the  movement  of  the  prototype  current  very  well. 

This  should  be  expected  since  it  was  entirely  consistent  with  the  theory 
presented  in  Appendix  A. 

44.  For  two  of  the  flood  events  that  were  simulated,  14  November 
1978  and  24  October  1976,  the  nose  of  the  density  currents  moved  slower 
in  the  model  than  in  the  prototype  throughout  most  of  the  tests  (Fig¬ 
ures  2  and  3).  This  is  probably  due  to  the  existence  of  a  viscous  flow 
regime.  Typically,  intrusion  speed  tends  to  decrease  as  viscous  forces 
start  dominating  current  movement.  Using  observed  data  from  the  14  No¬ 
vember  1978  event:  h  =  8.5  ft  ,  Ap  =  0.00072  g/cc  ,  p  =  1.0  g/cc  , 
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V  =  10  ft  /sec  ,  and  £  =  26,000  ft  (Figure  2),  Equation  A6  yielded  a 
value  of  123  which  is  much  greater  than  the  10.4  limit  for  nonviscous 
flow  (i.e.  the  prototype  is  well  into  the  nonviscous  flow  regime). 

From  test  6  of  this  event,  observed  model  data  are  h.^,  =  0 . 1  ft  ,  Ap.j, 

=  0.00035  g/cc  ,  p  =  1.0000  g/cc  ,  v  =  lO'^  ft^/sec  ,  and  £  =  7.33  ft  , 
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Figure  1.  Density  current  travel  distance  versus  time 
for  the  29  March  1979  flood  event 
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re  2.  Density  current  travel  distance  versus  time 
for  the  14  November  1978  flood  event 
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Figure  3.  Density  current  travel  distance  versus 
time  for  the  24  October  1976  flood  event 


Equation  A6  yielded  a  value  of  5  for  the  model,  indicating  that  the  model 
inflow  is  in  the  viscous  flow  regime.  Similar  results  were  obtained  for 
the  24  October  1976  event.  Since  the  density  current  in  the  model  was 
viscous-flow  dominated,  it  would  necessarily  move  slower  than  the  current 
in  the  prototype. 

45.  It  must  be  pointed  out  that  the  DeGray  distorted  physical 
model  was  constructed  prior  to  acquisition  of  the  prototype  data  for 
these  density  currents.  The  model  was  not  designed  to  simulate  these 
particular  events.  However,  analyses  of  the  results  from  each  event 
indicated  that  the  design  was  compatible  with  the  large  flood  of  March 
1979  but  not  with  smaller  events  of  November  1978  and  October  1976.  Be¬ 
cause  the  29  March  1979  flood  event  was  very  large  relative  to  the  inflow 
of  the  other  two  events,  the  flow  would  remain  nonviscous  in  the  model 
and  prototype  for  a  much  longer  distance.  Theoretically,  only  a  large 
flood  should  reproduce  well  in  the  model  since  the  distortion  of  this 
model  is  not  within  the  criteria  established  in  paragraph  28.  The  mini¬ 
mum  distortion  factor  for  this  model  based  on  Equation  25a  and  25b  should 
be  106.7;  the  actual  distortion  factor  is  20. 

46.  For  the  14  November  1978  flood  event  (Figure  2),  a  trend  for 
density  current  travel  time  was  observed  in  the  model  and  prototype: 

the  speed  decreased  as  the  current  progressed  through  the  lake.  For  the 
two  other  events,  similar  trends  were  observed  in  the  model,  although  it 
was  not  evident  from  limited  prototype  data.  The  velocity  decrease  is 
due  in  part  to  three-dimensional  effects.  From  the  work  of  Manins 
(1976),  for  two-dimensional  flow 

C  =  0.85  (qN)’^^  (32) 


where  the  variables  and  stratification  conditions  are  defined  in  Appen¬ 
dix  A.  For  two-dimensional  flow,  the  unit  discharge  q  is  the  total 
inflow  divided  by  the  width  of  the  current 


where 


Q  =  total  volumetric  inflow  rate 
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where  q(r)  =  unit  discharge  at  radius  r  .  From  Manins'  relationship 
and  Equation  35,  the  intrusion  velocity  would  also  be  inversely  related 
to  radial  distance.  Thus,  the  decrease  in  the  speed  of  the  density  cur¬ 
rent  in  the  model  and  prototype  is  at  least  qualitatively  consistent 
with  theory. 

47.  As  can  be  seen  in  Figures  2-3,  there  is  a  large  variability 
in  model  observations  on  presumably  "replicate"  tests  of  the  events. 

The  most  probable  cause  of  the  differences  from  test  to  test  is  the  dif¬ 
ficulty  in  creating  identical  stratifications  for  each  test  coupled  with 
model  construction  techniques.  The  model  is  made  with  "stepped"  sides 
(Figure  5)  to  permit  observation  of  current  movement.  The  elevations  of 
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Figure  5.  Stepped  model  construction  for 
convenient  observation 


neutral  bouyancy  for  the  inflow  currents  of  14  November  1978  and  29  Octo¬ 
ber  1976  were  about  el  367  ft  and  el  364  ft  NGVD,  respectively.  A  "step" 
in  the  model  occurred  at  el  370.  In  relation  to  current  thickness,  the 
differences  in  the  elevations  of  neutral  bouyancy  and  step  elevation  are 
very  small.  If  there  were  only  a  slight  error  in  inflow  density  that 
raised  or  lowered  the  neutral  bouyancy  elevation  in  the  model,  the  den¬ 
sity  current  could  experience  dramatic  differences  in  flow  width,  thereby 


causing  variation  from  test  to  test.  Additionally,  when  modeling  a  spe¬ 
cific  inflow  event,  it  is  necessary  to  not  only  reproduce  density  differ 
ences  but  also  the  density  gradient  of  the  reservoir  stratification. 
However,  variability  in  replicate  tests  is  not  evident  in  Figure  1.  The 
flood  event  on  29  March  1979  had  a  very  large  peak  inflow.  Because  of 
the  large  flow  rate,  the  current  was  very  thick  with  an  elevation  of  neu 
tral  bouyancy  of  354.  This  tended  to  eliminate  or  mask  the  effects  of 
small  differences  in  vertical  placement  of  the  density  current. 

48.  It  can  be  concluded  that  thickness  and  resulting  speed  of  the 
density  current  can  be  affected  by  viscous  stresses,  density  structure, 
and  model  geometry.  This  does  not  imply  that  distorted  physical  model 
analyses  cannot  be  effectively  employed.  Large  flows  may  be  simulated 
with  a  good  degree  of  confidence,  but  small  flows  must  be  evaluated 
regarding  the  forces  dominating  flow  movement  to  appropriately  design 
the  physical  model  to  be  used  in  their  study. 

49.  GRH  simulation.  Data  from  the  WES  Generalized  Reservoir 
Hydrodynamic  (GRH)  Facility  was  used  as  the  prototype  to  make  the  second 
model-prototype  comparison.  The  GRH  is  80.0  ft  long  and  2.85  ft  wide 
with  a  maximum  depth  of  3.0  ft.  Details  of  the  GRH  are  provided  in 
Figure  6. 


PLAN  VIEW 


50.  Two  conditions  were  simulated  in  the  GRH;  an  interflow  den¬ 
sity  current  and  an  overflow  density  current.  For  both  conditions,  in¬ 
flow  was  set  equal  to  outflow  with  a  flow  rate  of  10.0  gpm.  The  loca¬ 
tion  of  the  outflow  is  shown  for  the  two  conditions  in  Figure  6.  The 
inflow  was  baffled  to  distribute  the  flow  uniformly  from  top  to  bottom. 
Supply  water  for  the  GRH  can  be  heated  or  chilled  to  simulate  desired 
density  differences.  For  the  interflow  condition,  thermal  stratifica¬ 
tion  was  established  in  the  GRH  by  spreading  room  temperature  water  on 
top  of  chilled  water.  Room  temperature  water  was  used  for  the  upper 
layer  of  the  ambient  stratification  in  the  interflow  test  to  minimize 
surface  heat  transfer.  An  isothermal  condition  with  chilled  water  was 
set  up  for  the  overflow  case.  Temperature  profiles  are  shown  in  Fig¬ 
ure  7  for  the  two  conditions.  Chilled  water  was  warmed  and  automati¬ 
cally  controlled  to  provide  a  steady-state  inflow  temperature  of  62.0°  F 
for  the  interflow  condition.  Room  temperature  inflow  water  at  72.0°  F 
was  used  for  the  overflow  case.  This  minimized  the  surface  heat  trans¬ 
fer  of  the  overflow  density  current.  However,  since  the  ambient  water 
was  chilled,  it  was  possible  for  the  ambient  water  to  warm  slightly  dur¬ 
ing  the  test. 

51.  The  inflows  were  dyed  so  intrusion  length  and  thickness  (at 
various  stations)  could  be  visually  determined.  Additionally,  the 
elapsed  time  from  the  beginning  of  inflow  and  outflow  was  recorded.  The 
same  type  of  data  was  collected  from  the  distorted  model  so  that  a  direct 
comparison  could  be  made  with  the  prototype  (GRH). 

52.  The  prototype-to-model  horizontal  and  vertical  length  scale 
ratios  were  15.0  and  1.5,  respectively,  giving  a  distortion  factor  of 
10.0.  The  model  was  stratified  with  saline  and  fresh  water  to  provide 
tfie  same  density  difference  as  th.it  causeil  by  thermal  stratification  for 
the  GRH  interflow  rondition.  Simil.irly  the  salinity  of  the  model  inflow 
was  adjusted  to  provide  an  inflow  density  difference  (with  respect  to 
ambient  stratification)  erpia  I  to  Ih.it  of  the  GRH  interflow  condition. 

For  the  model  overflow  case,  a  freshwater  inflow  and  saline  receiving 
water  were  used,  again  maintaining  prototype  ilensity  difference. 

55.  The  two  inflow  conditions  were  simulated  with  the  model  and 
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two  types  of  comparisons  were  made.  Density  current  travel  time  or  in¬ 
trusion  length  versus  time  was  compared  for  the  model  and  GRH,  and  in¬ 
trusion  vertical  thickness  was  compared  at  various  longitudinal  stations 
for  the  same  intrusion  length.  The  comparison  of  model-prototype  travel 
time  for  the  interflow  case  is  shown  in  Figure  8.  Because  of  difficulty 
in  establishing  the  origin  for  elapsed  time  in  the  GRH  and  model  (due  to 
unequal  inflow  entrance  conditions),  all  times  are  referenced  to  the 
time  when  the  intrusion  reached  the  1/8  station  (or  traveled  10.0  ft  of 
the  GRH  length).  The  model  travel  time  compares  closely  with  that  of 
the  GRH  except  when  the  intrusion  reaches  the  end  of  the  flume.  It  is 
not  clear  why  the  GRH  interflow  slows  down  at  the  dam.  This  may  be  due 
either  to  the  beginning  of  viscous  effects  or  to  an  internal  wave 
(seiche)  reflecting  off  the  dam  (flume  wall)  as  the  interflow  is  arriv¬ 
ing,  thus  slowing  the  interflow  speed.  Seiche  action  in  the  model  may 
be  damped  or  out  of  phase  with  the  GRH.  In  general,  scale  models  with 
length  scale  distortion  can  be  used  to  model  long  waves  but  cannot  be 
used  to  model  short  waves  (LeMehaute  1976).  Thus,  it  is  possible  for 
internal  waves  of  short  length  to  be  out  of  phase  in  the  distorted  model 
compared  with  the  GRH.  The  GRH  and  model  intrusion  elevation  and  thick¬ 
ness  (Figure  9)  are  similar.  The  major  difference  between  them  is  that 
the  GRH  interflow  experienced  a  peak  thickness  where  the  model  interflow 
had  a  dip  (at  mid-length).  Again  the  difference  may  have  been  due  to 
unequal  damping  of  internal  waves  and  some  internal  wave  phase  shift. 

54.  Inspection  of  Figure  8  indicates  the  GRH  interflow  is  a  non- 

viscous  intrusion  for  £  <  70  ft  .  Solving  the  following  equation 

(Equation  A6)  for  h  0.275  ft  ,  Ap_,  =  0.00163  g/cc  ,  p  =  1.0  g/cc  , 
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H  =  70.0  ft  ,  and  v  =  10  ft  /sec  ,  results  in 


=  12.9 


(36a) 


Since  this  value  is  greater  than  10.4,  the  intrusion  is  in  the  nonviscous 
regime.  For  the  distorted  model 
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distance,  et 

Figure  9.  Interflow  density  current  thickness 

=71.1  (36b) 

for  the  same  intrusion  length.  Therefore  the  model  intrusion  is  also 
well  into  the  nonviscous  range.  This  is  not  surprising  since  the  dis¬ 
tortion  factor  was  greater  than  that  specified  by  Equation  25.  Clearly, 
the  intrusion  speed  and  thickness  were  adequately  reproduced  in  the 
model  as  both  the  model  and  GRH  intrusions  were  in  the  nonviscous  range 
for  most  (at  least  £  <  70.0  ft)  of  the  test. 

55.  It  is  worth  noting  that  the  use  of  Equation  M  (paragraph  A6 , 
Appendix  A)  to  predict  the  intrusion  thickness  resulted  in  h.j,  =  0.31  ft 
compared  to  h.^,  observed  of  about  0.275  ft.  Using  the  observed  Ap.j, 
of  0.00163  g/cc  and  h.j,  of  0.275  ft,  the  intrusion  speed  computed  from 
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Figure  11.  Overflow  density  current  thicknesses 


a  distance  of  60  ft,  the  GRH  travel  time  plot  curves  upwards  possibly 
due  to  viscous  effects.  The  curve  is  almost  straight  for  the  first  60  ft 
indicating  that  the  flow  must  be  in  the  nonviscous  range  for  £  <  60  ft. 
This  suggests  that  the  value  for  the  right  side  of  Equation  A9  may  be 
too  high.  The  model  overflow  was  computed  to  be  in  the  nonviscous  range 
throughout.  The  density  current  thickness  h  ,  as  computed  from  Equa- 
tion  A8,  is  0.14  ft  compared  with  about  0.15  ft  observed.  Use  of  Equa¬ 
tion  A7  results  in  a  computed  c  of  0.067  ft/sec  compared  with  that 
observed  of  c  =  40/470  =  0.085  ft/sec  .  This  fairly  good  agreement 
indicates  that  Equations  A8  and  A7  may  provide  a  good  first  approximation 
of  overflow  thickness  and  speed  for  the  nonviscous  range.  Again  the  dis¬ 
torted  model  agrees  relatively  well  with  the  GRH  since  the  physical  pro¬ 
cess  is  primarily  Froude  number  related  and  the  model  is  Freudian. 
Hydraulic  destratification 

58.  A  design  technique  for  lake  hydraulic  destratification  was 
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reported  by  Dortch  (1979)  and  Holland  and  Dortch  (1984).  In  hydraulic 
destratification,  water  jets  are  used  to  induce  mixing  (Figure  12). 


HYPOLIMNION 


f  f  y 


Figure  12.  Cross-sectional  view  of  a  reservoir  hydraulic 
destratification  system 

Although  there  are  several  physical  processes  that  contribute  to  mixing, 
such  as  entrainment  into  a  buoyant  jet/plume  and  subsequent  spreading  of 
the  flow  as  an  intrusion,  it  was  found  that  the  rate  of  mixing  could  be 
expressed  as  a  function  of  the  densimetric  Froude  number,  or 


' « 

^  =  0.20. 
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where 


t-  „  =  elapsed  pumping  time  to  reach  80  percent  totally 
destratified 

Q  =  pumping  rate 

V  =  total  lake  or  reservoir  volume 
V  =  exit  velocity  of  water  jets 
Ap  =  density  difference  of  hypolimnion  and  epilimnion 
p  =  density  of  epilimnion 
d,  =  depth  from  lake  water  surface  to  diffuser 


This  result  was  determined  through  an  experimental  program  at  WES  using 
laboratory  tanks  (see  paragraphs  71-72).  Turbulent  flow  was  maintained 
for  all  tests,  and  the  flow  rate,  jet  velocity,  and  reservoir  depth, 
volume,  shape,  and  density  stratification  were  varied  over  a  wide  range 
so  that  the  results  should  be  applicable  to  the  field  condition. 

59.  A  study  was  conducted  to  determine  if  the  effectiveness  of  a 
design  using  this  type  of  destratification  system  could  be  evaluated  on 
a  site-specific  basis  with  a  hydraulic  model.  If  so,  site-specific 
model  studies  to  improve  candidate  systems  would  be  desirable  prior  to 
field  installation.  Of  course,  for  practicality  and  preservation  of 
turbulent  flow,  the  lake  model  would  have  to  be  vertically  exaggerated. 

A  primary  concern  was  whether  the  mixing  process  could  be  modeled  pro¬ 
perly  with  a  distorted  model.  It  was  thought  that  a  vertically  exag¬ 
gerated  froudian  model  could  properly  simulate  this  type  of  destratifi¬ 
cation  because  the  mixing  process  was  dependent  on  a  densimetric  Froude 
number.  However,  the  mixing  process  depends  on  near-field  entrainment 
and  far-field  convective  spreading  so  there  was  some  question  as  to 
whether  a  distorted  model  could  properly  model  these  effects  simultan¬ 
eously  and  how  to  size  the  model  diffuser. 

60.  Because  there  were  no  field  data  on  this  type  of  mixing  pro¬ 
cess,  tests  previously  conducted  in  a  laboratory  tank  had  to  be  con¬ 
sidered  prototype  data  and  a  distorted-scale  model  of  the  tank  was  built 
Two  tests  conducted  with  a  20-f t-square ,  2.0-ft-deep  tank  were  selected 
for  study.  The  details  of  these  tests  are  summarized  below: 


whe  re  D 


Parameter 

Test  No.  1 

Test  No.  17 

Q,  gpm 

2.0 

3.0 

D ,  in . 

1/2 

1/2 

V,  fps 

3.27 

4.90 

Vl,  ft^ 

520.0 

720.0 

df,  ft 

1.3 

1.8 

Ap,  g/cc 

0.0036 

0.0022 

is  the  port 

diameter  of  the  single-port 

diffuser  that  was 

•  •  •  • 


44 


61.  A  distorted  model  of  the  tank  was  constructed  at  a  scale  of 


=  20.0  and  Y  =  2.0  .  As  before  Ap_  =  p  =  1.0  .  V  and  Q  were 
K  K  K  K  m  01 

scaled  according  to  Froudian  scaling  criteria  and  was  sized  to  give 

the  proper  V  given  Q  For  the  model,  the  properly  scaled  test 
m  m 

conditions  were: 


Parameter 

Test  No.  1 

Test  No.  17 

Q,  gph 

2.12 

3.18 

D,  in. 

5/64 

5/64 

V,  fps 

2.31 

3.46 

V,, 

0.65 

0.90 

h- 

0.65 

0.90 

Ap,  g/cc 

0.0036 

0.0022 

The  performance  of  the  mixing  systems  was  measured  by  plotting  the  per¬ 
cent  of  total  destratification  versus  pumping  time.  The  percent  of 
total  destratification  is  calculated  from  the  change  in  the  stability 
of  the  stratification  which  is  obtained  by  density  measurements  as  ex¬ 
plained  by  Dortch  (1979).  Because  the  percent  of  total  destratification 
is  dimensionless,  it  is  not  necessary  to  scale  those  numbers  from  model 
to  prototype.  Of  course,  model  pumping  time  had  to  be  scaled  to  proto¬ 
type  values. 

62.  The  performance  plot  for  test  No.  1  is  shown  in  Figure  13 
with  model  and  prototype  values  indicated.  It  is  apparent  from  this 
plot  that  the  model  simulated  the  prototype  very  closely.  When  model 
test  No.  17  was  run,  however,  it  was  found  that  the  model  mixed  much  too 
quickly.  The  model  depth  was  0.90  ft  for  this  test,  which  is  almost  the 
same  as  the  length  and  width  of  the  model  tank.  With  the  available  dis¬ 
tance  for  vertical  and  lateral  spreading,  it  was  considered  that  there 
might  be  wall  effects  that  could  influence  the  rate  of  mixing.  Thus,  a 
larger  vertical  scale  ratio  was  tried.  Using  the  model  depth  of  test 
No.  1  of  0.65  ft  for  test  No.  17  resulted  in  a  vertical  scale  factor  of 
Y,,  =  1.8/0.65  =  2.77  rather  than  Y„  =  2.0  .  This  resulted  in  a  DF 
of  7.22  rather  than  10.0.  The  No.  17  test  conditions  then  became: 


TIME.  MIN 


Figure  13.  Destratification  performance,  test  No.  1 


Parameter 

Value 

Q.  gph 

1.95 

D,  in. 

1/16 

V,  fps 

2.94 

v^, 

0.65 

dL,  ft 

0.65 

Ap,  g/cc 

0.0022 

This  adjustment  proved  to  be  successful  as  can  be  seen  by  the  close 
comparison  of  model  and  prototype  data  shown  in  Figure  14.  These  re¬ 
sults  indicate  that  this  type  of  destratification  system  can  be  simu¬ 
lated  fairly  well  with  a  Freudian  distorted  length  scale  model  as  long 
as  lateral  mixing  is  not  inhibited  by  the  model  walls.  If  the  distor¬ 
tion  factor  is  too  large,  the  model  walls  may  artificially  influence  the 
mixing  process. 
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PART  IV:  EXAMPLE  APPLICATIONS  OF  PHYSICAL  MODELS  OF  RESERVOIRS 

64.  Employing  much  of  the  information  discussed  in  previous  parts 
several  general  research  projects  and  site-specific  investigations  have 
been  conducted  at  WES  using  physical  models  of  reservoirs.  The  Bibli¬ 
ography  provides  a  publications  list  of  reports  on  studies  using  physi¬ 
cal  models  of  reservoirs.  Several  examples  are  described  below. 

Selective  Withdrawal  Research 

65.  Mathematical  descriptions  of  the  selective  withdrawal  phenome 
non  have  been  analytically  and  experimentally  developed  by  several  re¬ 
searchers  i.e.,  Bohan  and  Grace  (1969).  Simple  geometries  and  general 
(nonspecific)  conditions  have  been  and  are  being  used  to  develop  and 
enhance  the  existing  knowledge  about  selective  withdrawal.  Studies  on 
specific  project  models  have,  in  many  cases,  verified  or  enhanced  exist¬ 
ing  mathematical  destriptions  for  selective  withdrawal.  Research  de- 
signeil  to  refine  .uid  enlarge  the  applicability  of  existing  mathematical 
descriptions  continues  to  l>e  londucted  at  WES  using  small-scale  labora¬ 
tory  facilities.  I  he  it  I  mens  1  on  1  ess  relationships  developed  from  these 
studies  allow  the  .ipp  i  1 1  .i  t  i  on  of  model  results  to  prototype  situations. 

f)t) .  Selective  withdrawal  investigations  are  conducted  in  a  large 
flume  (Figure  IS)  designed  speiifically  for  selective  withdrawal  re¬ 
search.  St  ra  t  i  t  i  c  .1 1  1  on  similar  to  that  expected  in  reservoirs  and  lakes 
is  establisfied  fiy  slowly  flo.itiiig  fresh  water  over  saltwater.  This 
procedure  creates  tfie  epilimnion  and  hypolimnion.  Water  is  released 
th  rou  Rh  a  port,  over  a  weir,  or  via  other  outlet  devices  and  the  result¬ 
ing  velocity  distribution  and  withdrawal  zone  (Figure  16)  are  determined 
using  "dye  streak"  flow  visualization  techniques.  The  dye  streak  is 
displaced  toward  the  outlet  as  water  is  withdrawn  and  the  withdrawal 
<1  i st r i bu t i on  is  determined  from  the  time  rate  of  displacement.  The 
streak's  displacement  is  recorded  on  video  tape  allowing  stop-action 
analysis  of  the  withdrawal  zone  dimensions. 

67.  Ongoing  studies  are  being  conducted  to  evaluate  the 


•  _ 


Figure  16.  Selective  withdrawal  over  a  weir 


withdrawal  characteristics  of  various  outlet  configurations.  The  ef¬ 
fects  of  port  dimensions  (vertical  and  horizontal  extent  of  opening), 
p.irtially  submerged  ports,  and  submerged  sharp  crested  weirs  with  vari¬ 
ous  heads  and  weir  lengths  are  being  studied.  The  investigation  should 
result  in  the  development  of  criteria  that  define  the  flow  regime  over 
which  existing  selective  withdrawal  predictive  techniques  are  applicabl 
and  should  also  result  in  improvements  in  the  mathematical  descriptions 

Dickey-Lincoln  School  Lakes  Model 

68.  Dickey-Lincoln  School  Lakes  is  a  pumped-storage  hydropower 
project  proposed  by  the  New  England  Division  for  Maine.  An  undistorted 
near-field  model  of  Dickey  Lake  outlet  geometry  (outlet  structure  and 


local  topography)  was  constructed  to  study  the  withdrawal  and  pumpback 
mixing  characteristics  unique  to  the  Dickey  Lake  inlet-outlet  structure 
(Dortch,  Fontane,  and  Wilhelms  1976).  The  model  reproduced  the  structure 
and  upstream  topography  for  a  prototype  area  of  about  3,600  ft  square 
(Figure  17).  It  was  constructed  to  a  scale  of  1:200,  which  resulted  in 
the  following  scale  relationships  based  on  Froudian  modeling  criteria  for 
transferring  model  data  to  prototype  equivalents. 


Dimens  ions 

Scale 

Length  (prototype :model ) 

200 

Velocity 

14.14 

Time 

14.14 

Volume  flow  rate 

565,685 

Density  difference 

1 

69.  The  model  was  used  to  determine  if  the  vertical  distribution 
of  withdrawal  could  be  predicted  from  the  generalized  mathematical  se¬ 
lective  withdrawal  description  (Bohan  and  Grace  1973)  and  to  provide 
withdrawal  information  if  modifications  to  the  mathematical  description 
were  required.  The  model  was  also  used  to  evaluate  the  gross  entrain¬ 
ment,  mixing,  and  dilution  characteristics  of  the  pumpback  jet. 

70.  The  investigation  resulted  in  these  revisions  to  the  general¬ 
ized  mathematical  selective  withdrawal  description:  (a)  reduction  of 
effective  weir  length  to  force  a  lower  limit  of  withdrawal,  (b)  lowering 
of  the  predicted  elevation  of  maximum  velocity,  and  (c)  modification  of 
the  equation  that  predicts  velocity  distribution.  These  changes  were 
site  specific  and  attributed  to  the  effects  of  local  topography.  Pump¬ 
back  was  also  studied  in  the  model  and  entrainment  was  evaluated.  A 
computer  algorithm,  based  on  these  results,  was  developed  to  account 

for  the  mixing  characteristics  of  pumpback.  The  selective  withdrawal 
modifications  and  the  pumpback  algorithm  were  interfaced  with  the  one¬ 
dimensional  (vertical)  mathematical  thermal  reservoir  model  WESTEX 
(Holland,  Dortch,  and  Smith  1982)  for  simulation  of  the  spring  through 
fall  seasons. 
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Reservoir  Destratification 


71.  Artificial  destratification  of  reservoirs  is  one  of  several 
alternatives  for  increasing  the  dissolved  oxygen  (DO)  content  in  the 
lower  levels  of  a  reservoir.  Usually,  thermal/density  stratification 
in  a  reservoir  impedes  the  transport  of  oxygen  from  the  surface  to  the 
bottom  of  the  reservoir.  As  a  consequence,  the  hypolimnion  may  become 
anaerobic  resulting  in  dissolution  of  trace  metals,  release  of  nutrients 
that  may  stimulate  eutrophication,  formation  of  hydrogen  sulfide,  and 
depression  of  pH.  If  sufficient  mixing  energy  can  be  artificially  in¬ 
troduced  to  the  reservoir,  existing  stratification  can  be  eliminated 
allowing  complete  mixing  and  transport  of  DO  to  the  lower  levels  of  the 
reservoir,  thus  maintaining  DO  and  avoiding  the  effects  of  anaerobic 
conditions . 

72.  Research  has  been  conducted  by  Dortch  (1979)  and  Holland  and 
Dortch  (1984)  to  determine  the  design  parameters  and  provide  guidance  to 
produce  an  effective  means  for  artificially  destratifying  lakes  using 
hydraulic  methods  (pumping  water).  Laboratory  tanks  of  various  sizes 
and  shapes  were  used  to  simulate  stratified  reservoirs  to  develop  dimen¬ 
sionless  parameters  that  describe  the  destratification  process.  A  di¬ 
mensionless  description  of  the  mixing  phenomena  observed  in  the  model 
studies  was  developed  from  which  the  average  velocity  required  to  mix  a 
given  reservoir  in  a  specified  time  period  for  a  given  stratification 
and  pumping  rate  could  be  predicted.  These  results  can  be  used  in  the 
planning,  preliminary  design,  and  evaluation  of  a  hydraulic  system  for 
destratification  of  a  reservoir. 

Marysville  Lake  Hydrothermal  Study 

73.  This  study  was  conducted  to  determine  the  ability  of  the  pro¬ 
posed  Marysville  Lake  pumped-storage  project  (California)  to  satisfy 
downstream  water  temperature  objectives  (Fontane  et  al.  1977).  The 
study  required  (a)  an  understanding  and  description  of  the  large-scale 
hydrodynamic  phenomena  within  the  project,  (b)  development  of  expected 


temperature  profiles  within  Marysville  Lake,  and  (c)  estimation  of  tem¬ 
peratures  released  from  Marysville  Lake  and  from  the  afterbay.  A  physi¬ 
cal  hydrodynamic  model  and  a  numerical  simulation  model  were  used  in  the 
study.  The  distorted-scale  physical  model  provided  information  about  the 
hydrodynamic  response  of  Marysville  Lake  for  various  operational  condi¬ 
tions.  The  numerical  model  provided  the  capability  to  assess,  for  year¬ 
long  periods,  the  effect  of  historical  meteorologic  and  hydrologic  data 
and  various  operating  conditions  on  temperature  regimes  within  and  down¬ 
stream  of  the  project. 

74.  The  distorted-scale  model  was  used  to  evaluate  the  effects  of 
the  morphology  and  the  unsteady  generation  and  pumpback  operations  on 
the  density  stratification  of  Marysville  Lake  and  afterbay  and  to  help 
identify,  validate,  and  quantify  modifications  to  the  mathematical  model 
that  would  improve  the  reliability  of  the  predictions.  The  physical 
model  (Figurt^  18)  was  constructed  to  a  distorted  length  scale  ratio  of 
1:160  vertically  and  1:1,600  horizontally.  Through  vertical  scale  dis¬ 
tortion  it  was  possible  to  preserve  turbulent  flow  while  simulating  the 
entire  reservoir  and  afterbay.  Although  the  level  of  distortion  chosen 
did  not  conform  to  the  criteria  given  in  Equations  25a  and  b  for  the 
maintenance  of  the  nonviscous  spreading  regime,  the  processes  of  interest 
(i.e.,  entrainment)  were  dominated  by  inertial  rather  then  viscous  forces. 
Thus,  the  chosen  distortion  factor  did  not  compromise  the  model  results. 
With  the  density  differences  in  the  model  set  equal  to  those  in  the  pro¬ 
totype,  the  accepted  equations  of  hydraulic  similitude  based  on  Froudian 
relations  (paragraph  17)  were  used  to  express  the  mathematical  relations 
between  the  dimensionless  and  hydraulic  quantities  of  the  model  and  the 
prototype.  Allowing  for  vertical  scale  distortion,  the  general  relations 
for  transfer  of  model  data  to  prototype  equivalents  were  as  follows: 

_ Dimension _  Scale  Relation 

Length  in  vertical  direction  160 

(prototype/model ) 

Length  in  horizontal  direction  1,600 

(Continued) 


Dimension 


Area  in  vertical  plane 
Area  in  horizontal  plane 
Velocity 


Density  difference 


Scale  Relation 


256,000 

2,560,000 

12.65 

126.5 

3,238,172 


75.  It  was  determined  from  the  physical  model  study  that  the 
pumpback  density  current  had  a  much  stronger  effect  on  lake  stratifica¬ 
tion  than  inflow  or  withdrawal.  Inflow  to  the  reservoir  was  found  to 
affect  the  thermal  structure  of  the  lake  in  a  very  short  time,  usually 
within  a  week,  for  normal  river  flows.  It  was  also  determined  that,  due 
to  geomorphometry  of  the  reservoir,  a  large  volume  of  the  hypolimnion 
would  not  be  mixed  during  pumpback  and  would  not  be  released  during  gen¬ 
eration.  A  high  ridge  between  the  dam  and  the  Dry  Creek  arm  of  Marys¬ 
ville  Lake  prevented  the  cold  (hypolimnion)  water  in  the  Dry  Creek  arm 
from  being  mixed  or  withdrawn.  This  decrease  in  coldwater  storage  had 
to  be  accounted  for  in  the  numerical  simulations.  Additionally,  the 
amount  of  mixing  in  the  afterbay  was  estimated  with  the  physical  model 
and  these  results  were  incorporated  into  the  numerical  model.  The  modi¬ 
fications  made  to  the  numerical  model  helped  account  for  some  of  the 
hydrodynamic  and  geomorphological  effects  expected  to  occur  in  the  lake 
and  afterbay  and  provided  better  estimates  of  the  thermal  regimes  within 
and  downstream  of  the  project. 


PART  V;  SUMMARY 


76.  This  report  provides  guidance  for  the  use  of  physical  models 
to  study  reservoir  hydrodynamics.  Modeling  theory  is  developed  by  exam¬ 
ining  the  Navier-Stokes  equations  of  fluid  motion.  These  equations  are 
nondiraensiona 1 ized  resulting  in  fluid  motion  equations  containing  the 
Froude  and  Reynolds  numbers.  These  equations  are  used  to  develop  the  re¬ 
lationships  of  model-to-prototype  forces.  This  requires  that  the  model 
Froude  and  Reynolds  numbers  be  equal  to  prototype  Froude  and  Reynolds 
numbers,  respectively,  thereby  providing  complete  dynamic  similarity. 

77.  Upon  closer  scrutiny,  complete  dynamic  similarity  is  not 
feasible  if  ordinary  fluids  are  used  in  the  modeling  process.  In  other 
words,  if  water  is  the  model  and  prototype  fluid,  the  equivalence  of  the 
model  Froude  and  Reynolds  number  and  prototype  Froude  and  Reynolds  cannot 
be  preserved.  This  distortion  from  complete  similitude  does  not,  how¬ 
ever,  invalidate  physical  modeling.  Rather,  it  places  a  constraint  on 
modeling  that  must  be  considered  during  model  development.  For  free- 
surface  flow  regimes  such  as  reservoirs,  modeling  on  the  basis  of  the 
Froude  number  can  provide  accurate  similitude  if  the  model  Reynolds 
number  indicates  that  the  flow  is  of  the  same  fundamental  character 
(turbulent)  as  the  prototype.  Thus,  the  constraint  noted  in  the  previ¬ 
ous  paragraph  may  be  overcome  by  selecting  an  appropriate  length  scale 
for  the  model. 

78.  Scale  relationships  based  on  equating  the  model  and  prototype 
Froude  numbers  are  developed  for  undistorted-  and  distorted-length  scale 
models.  From  these  relationships,  the  ratio  of  prototype-quantity-to- 
model-quantity  can  be  determined.  Thus,  distance,  time,  velocity,  and 
flow  rate  may  be  scaled  for  transfer  from  model  to  prototype  or  proto¬ 
type  to  model. 

79.  Examples  of  applying  these  relationships  and  the  distorted- 
and  undistorted-scale  modeling  concepts  are  presented.  Several  studies 
are  briefly  discussed  and  important  results  of  the  physical  model  studies 
are  presented.  In  many  instances,  observed  hydrodynamic  processes  were 
mathematically  described  for  incorporation  into  a  numerical  hydrodynamic. 
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thermal,  or  withdrawal  model.  In  other  cases,  observations  improved  the 
credibility  of  numerical  predictions. 

80.  Undistorted-  and  distorted-length  scale  modeling  of  reservoir 
inflow,  outflow,  density  currents,  jets,  and  dispersion  is  discussed. 
Criteria  for  sizing  distorted  models  are  presented.  Effects  of  length 
scale  distortion  on  reservoir  hydrodynamic  processes  are  discussed. 
Recommendations  are  made  regarding  application  of  these  modeling  tech¬ 
niques  to  stratified  flow  phenomenon  in  reservoirs.  Froudian  models 
without  length  scale  distortion  should  be  used  to  physically  model  near¬ 
field  reservoir  hydrodynamics  which  would  include  studies  for  selective 
withdrawal  and  buoyant  jet  diffusion.  Physical  model  studies  of  far- 
field  reservoir  hydrodynamics  require  length  s,  lie  distortion  for  econ¬ 
omy  of  size,  preservation  of  turbulent  flow,  and  to  compensate  for  ex¬ 
cessive  interfacial  shear. 

81.  Viscous  and  inertial  forces  play  important  roles  in  the  move¬ 
ment  of  water  in  a  prototype  or  a  model.  Intuitively,  to  correctly  simu¬ 
late  these  forces,  a  model  should  reproduce  the  same  conditions  as  in  the 
prototype.  For  instance,  if  a  density  current  in  the  prototype  is  in  a 
nonviscous  flow  regime  (nonviscous  spreading),  then  the  model  current 
must  also  be  in  nonviscous  flow.  Criteria  were  developed  that  indicate 
whether  the  flow  is  viscous  or  nonviscous.  Distorted  model-prototype 
comparisons  are  made  for  density  current  advection.  The  results  veri¬ 
fied  that  nonviscous  spreading  of  a  density  current  in  the  prototype  can 
be  satisfactorily  simulated  with  a  Froudian  distorted-length  scale  model 
as  long  as  the  model  density  current  is  also  in  the  nonviscous  flow  re¬ 
gime.  It  is  also  demonstrated  that  hydraulic  destratification  of  a  res¬ 
ervoir,  which  is  a  Froudian  process,  can  be  modeled  very  satisfactorily 
with  a  distorted  scale  model. 

82.  In  many  cases,  results  from  near-field  flow  models  can  be 
coupled  with  far-field  distorted-length  scale  models;  often,  however, 
it  is  necessary  to  make  adjustments  in  the  distorted  model  to  the  near¬ 
field  process.  Furthermore,  when  modeling  far-field  flow  processes  that 
are  dominated  by  viscous  effects  and  turbulent  diffusion,  Froudian 
distorted- length  scaling  Is  necessary  but  will  not  ensure  similitude. 
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APPEtfl)IX  A:  MECHANICS  OF  INTERFLOW  AND  OVERFLOW  DENSITY  CURRENTS 


Interflow 

1.  In  his  work  on  interflows  (intrusions),  Manins  (1976)*  intro¬ 
duced  a  neutrally  buoyant  inflow  (evenly  distributed  laterally)  at  the 
mid-depth  of  a  linearly  stratified  laboratory  tank  (Figure  Al).  His 
experiments  were  in  the  nonviscous  range  where  there  is  a  balance  of 
pressure  and  inertia  forces.  In  this  range,  the  observed  intrusion 
length  £  was  linearly  proportional  to  time  t  .  Thus,  the  nose 
velocity  or  intrusion  speed  c  was  constant  for  a  given  set  of  condi¬ 
tions.  Through  dimensional  analysis,  the  intrusion  speed  reduced  to 

c  a  (qN)''^^ 

where 

q  =  volumetric  inflow  rate  per  unit  width 
N  =  buoyancy  frequency, 

=  inflow  density  and  density  of  ambient  fluid  at  mid-depth 
z  =  length  in  vertical  direction 


Figure  Al.  Sketch  of  intrusion  into  a  stratified 
fluid  (Manins  1976) 

Manins'  experimental  data  and  analytical  results  indicated  that  the  con 
slant  of  proportionality  in  the  above  relationship  was  about  0.85.  The 
intrusion  thickness  near  the  source  was  also  constant  during  each  test, 
although  the  thickness  does  taper  toward  the  nose  of  an  intrusion. 

See  References  at  the  end  of  the  main  text. 


Dimensional  analysis  yielded  for  the  half  thickness 


I 


a 
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Manins'  data  indicated  that  h  varied  nonlinearly  with  (q/N) 


was  not  the  case  for  c  versus  (qN)  .  For  Manins'  thickness  data, 


h/[(q/N)  ]  varied  between  about  1.7  and  1.0,  where  h  is  half  the 


intrusion  thickness  near  the  source. 


2.  Manins  analytically  showed  that  there  are  unsteady  flow  re¬ 
gions  due  to  columnar  wave  disturbances  which  are  able  to  propagate  up¬ 
stream  faster  than  the  relative  flow  velocity  c  .  There  are  M  such 
wave  fronts  defined  by 


M  =  integer  part 


(.-) 


where  H  is  half  the  channel  depth.  The  number  of  wave  fronts  tend  to 
affect  the  intrusion  speed  and  thickness.  The  larger  N  is,  the  larger 
M  is.  As  M  increases  (for  different  conditions)  the  internal  Froude 
number  of  the  intrusion  decreases.  Manins'  experimental  data  supported 
this  result  as  shown  in  Figure  A2  where  the  internal  Froude  number  F^ 
defined  by 


F  =  — 
r  Nh 


is  plotted  against  M  .  The  curves  are  upper  and  lower  bounds  from  an 
analytical  flow  model.  If  the  flow  is  completely  steady  throughout 
without  the  columnar  disturbances,  and  the  disturbance  of  the  ambient 
density  gradient  above  and  below  the  instrusion  is  negligible,  the  F 


would  be  unity;  but  as  M  increases,  F^  decreases  from  unity.  For 


Manins'  experimental  results,  F^  varied  from  about  1.0  to  0.40.  Noting 


F  is  dependent  upon  c  and  h  ,  it  is  evident  c  and  h  are  influ¬ 


enced  by  columnar  waves.  Manins  concluded  that  h  was  more  strongly 


influenced  by  columnar  waves  than  c  as  h  was  not  linearly  related  to 
1/2 
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(q/N)  (paragraph  1  of  this  appendix).  Thus,  although  the  intrusion 
speed  and  thickness  are  steady  in  the  nonviscous  regime,  internal  waves 
affect  the  intrusion  speed  and  thickness  (primarily  the  thickness)  with 


*•  •'  ••  "J 


•  •  •  *  •  .  “1  ^ 


-V 


• .  v'-.V 


C-'iL'dkj!) 


-V  . 

*•  *>■ 


* 


v-v-'.-.  ,--v-  -  ' 

WrfVA'.V.-v’.v' V. 

V'  -  -  -  - 


•  ••••• 


Figure  A2.  Internal  Froude  number  as  a  function 
of  the  number  of  forward  propagating  columnar 
modes  (Manins  1976) 


the  net  effect  of  reducing  F^  for  higher  M  .  Manins  did  not  give  any 
criteria  for  when  the  intrusion  is  in  the  viscous  or  nonviscous  regime. 

3.  Manins'  results  were  for  a  linear  stratification.  A  two-layer 
stratification  with  an  inflow  at  the  interface  is  shown  in  Figure  A3. 


Figure  A3.  Two-layer  stratification 


Following  from  Kao  (1977)  the  intrusion  speed  c  of  this  type  of  inter¬ 
flow  condition  can  be  described  for  the  portion  above  the  interface  by 


=  (2  ^  gh  -  pgh^) 


where 


p^(z)  =  density  gradient  within  the  density  current 
p  =  density  of  upper  layer 

h  =  vertical  distance  from  elevation  of  neutral  buoyancy  to 
upper  limit  of  interflow 


-1 


Assuming  that  the  density  gradient  within  the  density  current  is  linear, 
Equation  A1  reduces  to 


=  1.0 


(A2a) 


’‘•'-''if'' 


Therefore,  for  the  two-layer  stratification,  the  of  the  intrusion 

is  approximately  unity.  If  the  thickness  of  the  density  current  is  sym¬ 
metrical  about  the  interface. 
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where 


Ap.j,  =  the  total  density  differences  of  the  two-layer  stratification 
h.j,  =  the  total  thickness  of  the  intrusion 


Equation  A2a  can  be  written  as 


A*'  \ 


[IllCTtl 


4.  The  stratification  of  reservoirs  is  characteristic  of  the  two- 
layer  stratification  in  the  sense  that  there  exists  a  well-mixed  epilim- 
nion  and  hypolimnion.  These  two  layers  are  joined  by  a  metal imnion  that 
has  a  continuous,  often  linear,  stratification.  This  type  of  stratifi¬ 
cation  complicates  any  method  for  determining  M  ,  ,  c  ,  and  h  . 

As  a  first  approximation,  one  can  assume  that  reservoir  stratification 
is  more  characteristic  of  two  layers,  ignore  the  effect  of  internal 
waves,  and  assume  the  of  the  intrusion  is  unity.  Thus,  for  lake 

stratification.  Equation  A2  would  be  used  to  compute  c  if  h.p  is 
known.  Figure  A4  demonstrates  how  the  variables  of  Equation  A2  would  be 
applied  to  a  thermocline-type  stratification. 
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Figure  A4.  Definition  sketch  for  thermocline-type  stratification 

5.  Assuming  that  Equation  A2b  can  be  used  to  estimate  intrusion 
speed,  one  must  first  know  h.j,  to  use  the  equation.  The  thickness  of 
the  intrusion  is  nonuniform  and  decreases  from  h.j,  near  the  upstream 
source  to  zero  at  the  tip  of  the  nose.  With  q  constant,  the  average 
flow  velocity  increases  as  thickness  decreases  with  the  maximum  veloc¬ 
ity  c  existing  at  the  nose.  Since  the  flow  is  nonuniform,  the  use  of 
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q  =  ch^  is  not  the  appropriate  relationship  to  couple  with  Equation  A2 
to  solve  for  c  and  h^  .  However,  it  is  true  that  q  is  proportional 
to  ch^  and,  thus,  after  rearranging  and  making  use  of  Equation  A2, 


(A3) 


Relationship  A3  can  also  be  obtained  from  h^  «  (q/N) 


1/2 


(paragraph  1) 


by  approximating  N  as  ■^g/p(Ap^/h^)  .  A  study  of  Manins'  data  indi¬ 
cated  that  the  proportionality  constant  for  relationship  A3  varied  be¬ 


tween 


:.5  I  h/  [(q/N)^''^J  =  l.oj 


]  =  ■ 


and  about  5.1  ^h/  j(r 

with  the  variability  due  to  variations  in  M  (paragraphs  1  and  2).  Us¬ 
ing  approximately  the  mean  value  of  Manins'  data  for  the  thickness  rela¬ 
tionship  results  in 


2/3 


=  3.0 


(A4) 


for  a  symmetrical  interflow.  With  Ap.j.  dependent  upon  h.j,  ,  h.p  can 
be  estimated  from  an  iterative  solution  of  Equation  A4.  With  a  value 
for  h.j.  ,  c  can  be  estimated  (for  reservoir  stratification)  from  Equa¬ 
tion  A2.  The  reservoir  width  used  for  computing  q  (for  Equation  A4) 
would  be  taken  near  the  station  at  which  the  inflow  separates  from  the 
bottom  and  becomes  an  interflow. 

6.  The  results  of  Manins  can  be  used  for  an  intrusion  that  is 
asymmetric  about  the  neutrally  buoyant  layer  as  follows: 

2/3 


^1,2  = 


(A5) 


fAp 
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hj  2  ~  the  intrusion  thickness  below  (1)  and  above  (2)  the 
’  neutrally  buoyant  layer 


where 


the  density  difference  between  the  elevation  of  lower 
limit  (1)  or  upper  limit  (2)  of  the  intrusion  thickness 
and  the  elevation  of  the  neutrally  buoyant  layer 


It  should  be  emphasized  that  for  the  the rmoc line -type  stratification 
the  will  be  less  than  1.0  and  therefore  the  use  of  Equation  A2  for 

interflow  speed  will  provide  only  an  estimate.  If  a  linear  stratifica¬ 
tion  was  observed,  Manins'  results  could  be  used  to  obtain  a  more  accu- 

1/2 

rate  estimate.  For  the  linear  case,  c  =  0.85(qN)  ,  where  N  is 

known  and  h  is  not  required  to  find  c  .  After  calculating  c  ,  the 
values  for  M  ,  determined  (paragraph  2).  But  for 

the  arbitrary  thermocline-type  stratification,  N  is  not  known,  and  c  , 
M  ,  F^  ,  and  h  cannot  be  calculated  directly  from  Manins'  results; 
thus,  the  approximations  above  are  proposed.  It  must  also  be  emphasized 
that  the  above  discussions  pertain  only  to  two-dimensional  intrusions. 

7.  Imberger,  Thompson,  and  Fondry  (1976)  suggest  that  an  interflow 
intrusion  is  in  the  nonviscous  range  when 


t  <  su 

-  Nv 


where  t  is  the  elapsed  time  from  when  the  intrusion  began.  Thus 
£  =  ct  is  the  intrusion  length.  Solving  Equation  A4  for  q  ,  using 
t  =  £/c  and  N  =  and  substituting  Equation  A2b  for  c 

Imberger' s  relation  can  be  rewritten  as 


>  10.4 


In  this  form  Equation  A6  represents  the  product  of  a  hydraulic 
grade  line  h_/£  and  the  densimetric  Froude-Reynolds  number 


l(^^/p)g/\j 


Equation  A6  provides  criteria  for  cases  where 


a  nonviscous  interflow  exists. 


Overflow 


8.  Koh  (1976)  summarizes  the  results  for  a  variety  of  buoyancy- 
driven  gravitational  surface  flows.  His  results  for  the  two-dimensional 
continuous  surface  inflow  pertain  mostly  to  reservoir  overflow  density 
currents.  For  this  case,  Koh  presented  analytical  and  experimental  re¬ 
sults  that  compared  closely  and  related  the  speed  of  advance,  for  inter¬ 
mediate  travel  times  (nonviscous  regime),  as 


*■>) 


1/3 


where  Ap  is  now  the  density  difference  of  the  inflowing  water  and  the 
receiving  water.  The  above  equation  reduces  to 


c 


1/3 


gq 


(A7) 


9.  From  Koh's  (1976)  results  it  was  concluded  for  a  surface  cur¬ 
rent  that 


2/3 


(A8) 


g 


where  h^^  is  the  thickness  of  the  overflow  density  current.  Equation  A8 
is  fairly  consistent  with  Equation  A5 ,  the  upper  or  lower  thickness  equa¬ 
tion  for  an  interflow.  If  Equation  A8  is  solved  for  q  and  substituted 
for  q  in  Equation  A7,  Equation  A7  can  be  rewritten 


c  =  0.90  ^^h^ 


which  is  close  to  the  previous  interflow  result  of  F^  near  unity  (para¬ 
graphs  2  and  3). 

10.  It  was  observed  from  Koh's  (1976)  results  that  the  surface 
intrusion  length  was  linearly  proportional  (nonviscous  regime)  to  time 
for  less  than  approximately  0.5  to  1.0.  For  greater 


A8 


than  this  range,  the  intrusion  length  was  proportional  to  t  (viscous 
regime).  Thus,  the  criterion  for  a  nonviscous  overflow  is 

^  <  0.5  to  1.0 
o 


where  the  parameter  t^  is  defined  as 


t 

o 


16q^ 

3n^v 


Substitution 
and  using  t 
results  in 


of  this  relation  for  t^  into  the  previous  relationship, 
=  £/c  with  Equation  A7  for  c  and  Equation  A8  for  q 


t  ^ 


HV 


>  4.88  to  9.75 


(A9) 


The  larger  value  for  the  right  side  of  Equation  A9  is  very  close  to  the 
value  in  Equation  A6.  A  range  of  values  for  the  right  side  of  Equa¬ 
tion  A9  is  presented  because  the  data  reported  by  Koh  (1976)  gradually 
transition  from  the  nonviscous  to  the  viscous  regime. 

11.  To  summarize,  the  following  equations  may  be  used  for  esti¬ 
mating  interflow  and  overflow  conditions: 


Interflows 


Equation 


Nonviscous  range 

Thickness 

Speed 


A6 

A4  and  A5 
A2 


Overflows  Equation 


Nonviscous  range  A9 
Thickness  A8 
Speed  A7 


In  the  nonviscous  regime,  the  thickness  near  the  source  and  the  speed  of 
the  intrusion  are  constant.  In  the  viscous  range,  the  intrusion  speed 


decreases  with  time  and  the  thickness  increases  with  time.  The  viscous 
regime  will  not  be  discussed  in  any  greater  detail  because  it  is  not 
practical,  and  most  difficult,  to  physically  model  with  Froudian  scaling 
criteria  (see  Part  III  of  the  main  text). 


.ns: 


APPENDIX  B:  PROCEDURE  FOR  SIZING  A  ROUND  BUOYANT  JET 
IN  A  DISTORTED  MODEL 


1.  If  the  entrainment  (dilution)  of  a  buoyant  jet  is  known,  it  is 
possible  to  size  the  jet  in  the  distorted-scale  physical  model  to  give 
the  same  dilution.  A  round  buoyant  jet  discharging  horizontally  into  a 
linearly  stratified  reservoir  has  an  ultimate  (finite)  dilution  that 
depends  upon  the  specific  physical  conditions.  Roberts  (1982)*  reported 
procedures  for  determining  the  dilution  for  this  type  of  jet.  He  ex¬ 
panded  upon  the  results  of  Fan  and  Brooks  (1969)  to  include  the  condi¬ 
tions  of  a  weakly  buoyant  jet.  These  results,  which  include  weakly  and 
strongly  buoyant  jets,  are  discussed  and  the  procedure  for  properly  siz¬ 
ing  (scaling)  a  round  jet  is  presented  in  subsequent  paragraphs. 

2.  Consider  the  following  known  parameters  for  a  prototype  jet: 
the  initial  (source)  volume  flux  of  the  jet  before  entrainment  Q  ,  and 
the  area  of  efflux  A  ,  the  difference  in  density  between  the  jet  and 
ambient  Ap^  ,  and  the  stratification  parameter  s  .  Mathematically 
expressed 


Ap  =  p  -  p 

o  o 


e  =  -g/p  dp/dz 


where 


p^  =  initial  density  of  jet 
p  =  ambient  density  at  the  level  of  jet  efflux 
g  =  acceleration  due  to  gravity 
dp/dz  =  slope  of  ambient  density  gradient  at  level  of  efflux 
With  these  prototype  values,  it  is  possible  to  calculate  the  exit  ve¬ 


locity  u^  ,  momentum  flux  M  ,  and  buoyancy  flux  B  as  follows: 


o  A 


M  =  Qu 


See  References  at  the  end  of  the  main  text 
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3.  By  combining  £  ,  Q  ,  M  ,  and  B  ,  characteristic  lengths 
can  be  formed  that  describe  jet  properties.  The  most  convenient  are: 


0  =  - 


..  •  © 


where  is  proportional  to  the  distance  over  which  the  nozzle  geom¬ 
etry  is  important  in  determining  the  jet  properties;  is  proportional 

to  the  distance  over  which  the  initial  momentum  flux  is  important  rela¬ 
tive  to  the  initial  buoyancy  flux;  and  is  proportional  to  the  dis¬ 

tance  over  which  the  source  momentum  flux  is  important  relative  to  the 
stratification.  Jet  properties  and  characteristics  can  be  expressed  as 
ratios  of  these  lengths. 

4.  Through  dimensional  analysis,  total  volume  flux  p  in  the  jet 

following  entrainment  can  be  related  to  and  in  the  form 

Q  e  M  e 

1/4  [£-  £„ 

-  =  f  -2  — y  ( B 1  ) 

3/4  ^  £  ’  £ 

e  £ 


The  dimensionless  parameters  of  Equation  B1  are  plotted  in  Figure  B1  for 
the  case  of  a  round  buoyant  jet  discharging  horizontally  into  a  linearly 
strafified  ambient  as  reported  by  Roberts  (1982). 

5.  From  Figure  B1  it  is  possible  to  determine  p  for  the  proto¬ 
type  jet,  and  the  entrainment  parameter  E  which  is  related  to  p  by 

_  p  -  Q  _  ent ra i ned  volume  flux 
Q  initial  volume  flux 


Mean  dilution  of  the  jet  is  defined  as  p/Q  . 
of  the  prototype  jet  in  the  model  requires 


1.0 


To  preserve  the  dilution 
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Figure  Bl.  Dimensionless  volume  flux  parameter  for  round  buoyant 
jets  discharging  horizontally  into  a  linearly  stratified  ambient 


or 


where  R  represents  the  prototype-to-model  ratio  of  the  parameter.  With 

p  (p  =  prototype),  p  (m  =  model)  can  be  calculated  from  Equation  B2 . 

pm 

Since  the  discharge  ratio  Q  is  scaled  according  to  Froudian  criteria 
3/2 

^R  ~  ^R^R^  '  Therefore, 


M 


m 


X 

^r’'r 


(B3) 


where  X.  ,  Y  =  prototype-to-model  ratio  for  length  in  horizontal  and 
H  R 

vertical  directions,  respectively.  If  the  water  density  differences  of 
the  model  stratification  are  equal  to  that  of  the  prototype,  then  the 
stratification  parameter  is  scaled  as 


•  •  •  •  • 


•  •  •  •  • 


B3 


The  momentum  flux  for  the  model  can  then  be  written 


2  2 

M  n 

M  =  Q  u  =  ^  =  o — 
m  mo  A  ^2. 

m  m  Qr,A 
R  P 


where  =  area  scale  ratio  in  a  vertical  plane.  Combining  Equa¬ 

tions  B3  through  B5 


m  \  /  p 


1/2 

R 


where  is  the  prototype-to-model  ratio  of  the  jet  source  diameter. 

R  r  •%  1 1 


The  parameter  was  determined  previously  from  Fig- 

*■  ^  r  1/A  3/A  I 

ure  Bl.  Equation  B6  has  two  unknowns  (|je  )/(M  )J  ^  and  . 

6.  The  Jet  characteristic  length  ratios  of  paragraph  3  of  this 


appendix  can  be  expanded  to  the  form 


0  3/A  1/4 

_2  =  A  £ 
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From  the  above  relations  with  Ap_  -  1.0  ,  the  following  relationships 

K 

result 
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